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ABSTRACT 

Extensive analyses of statistical strong gravitational lensing are performed based 
on the final Cosmic Lens All-Sky Survey (CLASS) well-defined statistical sample 
of flat-spectrum radio sources and current estimates of galaxy luminosity functions 
per morphological type. The analyses are done under the assumption that galactic 
lenses are well-approximated by singular isothermal ellipsoids and early-type galaxies 
evolved passively since redshift z ^ 1. Two goals of the analyses are: (1) to con¬ 
strain cosmological parameters independently of other techniques (e.g. Type la super¬ 
novae magnitude-redshift relation, cosmic microwave background anisotropies, galaxy 
matter power spectra); and (2) to constrain the characteristic line-of-sight velocity 
dispersion and the mean projected mass ellipticity for the early-type galaxy popu¬ 
lation. Depending on how the late-type galaxy population is treated (i.e., whether 
its characteristic velocity dispersion is constrained or not), we find for a flat uni¬ 
verse with a classical cosmological constant that the matter fraction of the present 
critical density Dm = 0.3lto)i4 (68%) for the unconstrained case or 0.40to)i6 (68%) 
for the constrained case, with an additional systematic uncertainty of « 0.11 arising 
from the present uncertainty in the distribution of CLASS sources in redshift and 
flux density. For a flat universe with a constant equation of state for dark energy 
w = Pa; (pressure)/pa; (energy density) and the prior constraint w > —1, we find that 
— 1 < w < — 0 . 55 ^ 0)11 (68%) for the unconstrained case or —1 < w < —0.4llg ^g 
(68%) for the constrained case, where w = —1 corresponds to a classical cosmologi¬ 
cal constant. The determined value of the early-type characteristic velocity dispersion 
(cr**^) depends on the faint-end slope of the early-type luminosity function and 

the intrinsic shape distribution of galaxies; for equal frequencies of oblates and pro- 
lates, we find that ct1'^^(0.3 < 0 < 1) = IQStig (68%) for a ‘steep’ = —1 or 

(0.3 < z < 1) = 1811]^5 km s“^ (68%) for a ‘shallow’ = —0.54. Finally, from 
the relative frequencies of doubly-imaged sources and quadruply-imaged sources, we 
find that a mean projected mass ellipticity of early-type galaxies Cmass = 0.42 with a 
68% lower limit of 0.28 assuming equal frequencies of oblates and prelates. 

Key words: gravitational lensing - cosmological parameters - galaxies: structure - 
galaxies: kinematics and dynamics - galaxies: halos - methods: statistical 


1 INTRODUCTION 

Analysing the statistical properties of strong gravitational 
lensing in a sample of cosmologically distant sources pro¬ 
vides a method for constraining cosmological parameters 
(e.g. Turner 1990; Fukugita, Futamase, & Kasai 1990; 
Fukugita & Turner 1991; Fukugita et al. 1992; Carroll, Press, 
& Turner 1992) as well as for constraining global proper¬ 


ties of galaxy populations including evolutions (e.g. Maoz & 
Rix 1993; Kochanek 1993; Mao & Kochanek 1994; King & 
Browne 1996). The observed statistical properties of gravi¬ 
tational lensing in a sample are the total rate of multiple¬ 
imaging, the image separations, the lens redshifts, the source 
redshifts, and the image multiplicities. These properties de¬ 
pend not only on cosmological parameters but also on the 
following properties of galaxies and sources (e.g. Turner, Os- 
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triker, & Gott 1984; Turner 1990; Fukugita & Turner 1991; 
Mao 1991; Kochanek 1993; Maoz & Rix 1993; Wallington 
& Narayan 1993; Mao & Kochanek 1994; King & Browne 
1996; Kochanek 1996a,b): (1) the galaxy luminosity func¬ 
tions per morphological type (and their evolutions of any 
relevance to lensing) and the projected mass distribution in 
the inner cylindrical region of a galaxy of given luminos¬ 
ity and morphological type; and (2) the distribution of the 
sources in redshift and flux in the observational selection 
waveband (e.g. radio). Because of the dependences of the 
statistical properties of gravitational lensing on the above 
multiple factors, analyses of statistical lensing can be used 
for various purposes. On the other hand, in order to derive 
reliable results from an analysis of statistical lensing, it is 
vital to properly take into account all the factors through a 
statistical lensing model based on reliable data. 

In this paper we perform extensive analyses of statistical 
lensing based on more reliable sets of data than were avail¬ 
able previously and using a more realistic statistical lensing 
model. Essential results on cosmological parameters from 
(part of) the analyses are reported in Chae et al. (2002). 
However, in this paper we present not only the details of the 
analyses but also comprehensive results of the analyses for 
a broader range of parameters. In particular, we derive and 
discuss the global properties of the galaxies. In comparison 
to previous works on statistical lensing that attempted to 
put limits on cosmological parameters (e.g. Kochanek 1996a; 
Falco, Kochanek, & Munoz 1998; Chiba & Yoshii 1999; Hel- 
big et al. 1999; Cooray 1999), the following aspects of this 
work can be pointed out: 

(i) This work is based on the final Cosmic Lens All-Sky 
Survey (CLASS: Myers et al. 2003; Browne et al. 2003) sta¬ 
tistical sample of 8958 radio sources [which includes Jodrell- 
Bank-VLA Astrometric Survey (JVAS) sources as a sub¬ 
sample] containing 13 multiply-imaged sources, which is the 
largest statistical sample that satisfies well-dehned observa¬ 
tional selection criteria. 

(ii) This work is based on estimates of luminosity func¬ 
tions of galaxies per morphological type that have been de¬ 
rived taking into account latest results from large-scale ob¬ 
servations of galaxies particularly including the Two Degree 
Field Galaxy Redshift Survey (2dFGRS) and the Sloan Dig¬ 
ital Sky Survey (SDSS). 

(iii) We adopt a statistical lensing model in which we not 
only use elliptical projected densities for galaxies but also in¬ 
corporate intrinsic shapes of galaxies. The adopted galactic 
model is the axisymmetric (i.e. oblate and prolate) singu¬ 
lar isothermal ellipsoid (SIE) with a distribution function of 
the form T — T(E,Lz), where E and Lz are respectively 
the relative energy and the angular momentum component 
parallel to the symmetry axis. 

(iv) This work is intended to provide robust constraints 
in cosmological parameter space. We make the least pos¬ 
sible prior assumptions on the parameters of the assumed 
statistical lensing model and constrain the parameters self- 
consistently from the data of statistical lensing. 

(v) This work is intended to be an investigation of global 
properties of galaxy populations; namely, characteristic ve¬ 
locity dispersion(s) for the early-type (and the late-type) 
galaxy population(s), mean apparent mass axial ratio of 
galaxies, and intrinsic shapes of galaxies. These investiga¬ 


tions are made possible by the use of the (3-dimensional) 
ellipsoidal lens model [the above point (iii)]. 

This paper is organized as follows. In section 2, we 
build up the statistical lensing model based on the singular 
isothermal ellipsoid (SIE) mass model for a galaxy, includ¬ 
ing the calculation of the cross sections for double-imaging 
and quadruple-imaging as functions of line-of-sight velocity 
dispersion, apparent ellipticity and assumed intrinsic shape 
distributions, and the calculation of magnihcation biases for 
double-imaging and quadruple-imaging. In section 2, we also 
define the likelihood function for the data given the model of 
statistical lensing. In section 3, we summarize the input data 
for the analyses including the final CLASS statistical sam¬ 
ple, distribution of the CLASS sources in redshift and radio 
flux density, and galaxy luminosity functions per morpho¬ 
logical type. The results of fitting the model to the data are 
presented in section 4. In section 5, we discuss the derived 
results including comparison with previous results from sta¬ 
tistical lensing and examination of possible systematics that 
could affect the derived results. In section 6, we discuss fu¬ 
ture prospects of statistical lensing. Conclusions are given 
in section 7. 


2 MODEL AND METHOD 

2.1 Lens model: the singular isothermal ellipsoid 

The isothermal mass profile (i.e. the density p ~ ^~^) 
has been widely used for analyses of statistical lensing 
(e.g. Turner, Ostriker, & Gott 1984; Fukugita & Turner 
1991; Kochanek 1993; Maoz & Rix 1993; Wallington & 
Narayan 1993; Mao & Kochanek 1994; King & Browne 1996; 
Kochanek 1996a,b; Falco, Kochanek, & Munoz 1998; Chiba 
& Yoshii 1999; Helbig et al. 1999; Cooray 1999; Rusin & 
Tegmark 2001), because of the analytical tractability of the 
isothermal profile and because the projected surface density 
of the inner cylindrical regions of lensing galaxies probed by 
strong gravitational lensing are either consistent with or not 
too different from the isothermal profile.^ Isothermal mass 
models would, in general, be too simplihed to allow us to do 
accurate modelling of individual lenses (see, e.g., Munoz, 
Kochanek, & Keeton 2001; Chae 2002), but they should 
suffice at least as hrst-order approximations to the mean 
properties of galaxies relevant to statistical lensing. Thus, 
following the literature on statistical lensing, we adopt the 
isothermal mass proHle for our analyses. However, we use a 
more general isothermal model that incorporates not only a 
galaxy’s apparent axial ratio but also its intrinsic shape.^ 


^ The reader is referred to numerous examples in the literature 
[see, e.g., Treu & Koopmans (2002) and Koopmans & Treu (2003) 
for very recent examples]. 

^ All the previous works on statistical lensing that intended to 
put limits on cosmological parameters used an isothermal sphere, 
while the works that focused on studying the relative frequencies 
of image multiplicities used an elliptical lens model (e.g. King & 
Browne 1996; Kochanek 1996b; Rusin & Tegmark 2001). 
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2.1.1 The surface density, dynamical normalisations, and 
cross sections 


The projected surface density of the SIE lens S®™(a;,y) can 
be written in units of the critical surface density Scr by 

^cr 

D{0,zs) 

_ re. VfHf) 

2 y^a;2 + Py^ ’ 


where / is the minor-to-major axial ratio of the surface den¬ 
sity and is equal to 1 — e (where e is the ellipticity), Rh is 
the Hubble length cjHo, D{z\,Z 2 ) is the angular diameter 
distance between redshift zi and 22 (22 > 21 ) in units of 
Rh, and zi and Zs are respectively the lens and the source 
redshifts. In equation (1), r.. is the critical radius for the 
singular isothermal sphere (SIS) (i.e. for f ~ 1) and is given 
by 


r 


cr 


47r 



D(0,zi)D{zi,Zs) 

- 6 ( 0 , Zs) 


( 2 ) 


where a is the line-of-sight velocity dispersion of the galaxy.® 
For the case of the SIS, only double imaging is allowed 
and the cross section for the double imaging is Trrf. For 
/ 7 ^ 1 , quadruple imaging (and triple imaging for a suffi¬ 
ciently flattened surface density; see below) is possible in 
addition to double imaging, and the relations between the 
multiple-imaging cross sections and the line-of-sight veloc¬ 
ity dispersion of the galaxy are not unique but depend on 
the intrinsic shape, the viewing direction, and the distribu¬ 
tion function of the galaxy. This dependence is parametrised 
by A(/) (‘dynamical normalisation’) in equation (1). In our 
parametrisation of the SIF in equation (1), A(/ = 1) = 1 and 
the case A(/) = 1 for any / corresponds to the normalisa¬ 
tion in which the mass within the area enclosed by an equi- 
density ellipse is independent of / (see Kormann, Schneider, 
& Bartelmann 1994). As shown in Kormann, Schneider, & 
Bartelmann (1994), for / < /o(~ 0.3942) triple imaging 
with comparable brightnesses for all three images is possi¬ 
ble due to a naked cusp radial caustic. The cross section for 
the naked-cusp triple-imaging is expected to be small and no 
case has been found from the CLASS (Browne et al. 2003). 
Hence we ignore the possibility in the following analyses. 

Let S 2 {f) and S 4 (/) denote respectively the double 
imaging and the quadruple imaging cross sections in units 
of for the case of A(/) = 1 for any /. Then, for any 
given A(/), the double (quadruple) imaging cross section 
S 2 (/) [S4(/)1 in units of r'^, is given by S 2 {f) = [A(/)1^S2(/) 
= ['^(/) 1 ^'S 4 (/)}- The cross sections for the double 
imaging and the quadruple imaging for A(/) = 1 are respec¬ 
tively given by (Kormann et al. 1994) 


«4(/) 


4 / 

l _/2 


^ / Vl — a;® 


— arccos x 


4^4., (3) 


^ Here we are considering a single-component mass model for the 
total mass distribution of the galaxy. 


and 


Mf) = ar - 2 s4(/) 


( 4 ) 


which are valid for / > /o. Notice that the cross sections 
given by equations (3) and (4) are for all possible image 
magnification ratios of the multiple images. However, for the 
CLASS statistical sample (section 3.1) the observational se¬ 
lection limit on the image magnification ratio for the double 
imaging is 7?. > 0.1 where TZ is the ratio of the flux densities 
of the fainter to the brighter component. Thus, the double 
imaging cross section that satisfies the observational limit 
on the image magnification ratio is less than that given by 
equation (4) but it can be calculated numerically through 
a Monte Carlo method. The total cross section for multiple 
imaging for the SIE is given, in units of r.., by 


Stotif) = S2{f)+S4{f) 

= [A(/)]^stot(/) 

= [A(/)f[S2(/)+S4(/)]. 


( 5 ) 


The rest of this subsection is devoted to the calculation 
of the dynamical normalisation factor A(/) in equation (1). 
For the calculation of the dynamical normalisation, we need 
to calculate the average of the line-of-sight velocity disper¬ 
sions for an ensemble of models that are allowed by the axial 
ratio / of the surface mass density (equation 1). In gen¬ 
eral, calculations of velocity dispersions of a galaxy require 
numerically constructing the orbits of the particles in the 
galaxy that satisfy physical distribution functions IF > 0.^^ 
However, in this paper we consider the cases where the 
velocity dispersions can be calculated analytically. These 
are the axisymmetric models, namely the oblate and the 
prolate spheroid, with distribution functions of the form 
T = T(E, Lz), where E and Lz are respectively the relative 
energy and the angular momentum component parallel to 
the symmetry axis. The axisymmetric singular isothermal 
mass model can be written in the cylindrical coordinates 
[R, Z) and in the polar coordinates (r, P) by 


P = Po 




1 


o-Q 


+ g-2^2 27rG + g-2^2 
1 


27rG r2(sin® i9 + q~^ cos^ i7) 


( 6 ) 


where q is the intrinsic axial ratio of the spheroid. The oblate 
and the prolate spheroids respectively correspond to g < 1 
and q > 1. For q = 1 (i.e. the SIS case), the parameter ao is 
identical to the one-dimensional velocity dispersion a (equa¬ 
tion 2). Under the assumption that the distribution function 
has the form T = T(E,LP), the Z-component g\ and the 
-R-component a\ of the second velocity moments are iden¬ 
tical [see, e.g., Binney & Tremaine 1987, equation (4-175)]. 
The azimuthal component of the velocity distribution at a 
point includes streaming rotation motion as well as random 
motions. We use the mean squared azimuthal velocity, to 
be denoted by a'f, for our dynamical normalisation calcu¬ 
lations since measured line-of-sight velocity dispersions over 
a large aperture include contributions from any streaming 


® In this paper, we use the notation T for the distribution func¬ 
tion, since the notation / is used to denote the surface density 
axial ratio. 
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motions. For the axisymmetric scale-free isothermal model 
given by equation (6), the second velocity moments can be 
calculated analytically and the expressions include only ele¬ 
mentary functions (see, e.g., van der Marel 1994; Qian et al. 
1995). The second velocity moments are given for the oblate 
and the prolate cases as follows: (1) For g < 1, we have 


= u|(d) 

= cro(l + cot^'ll) 

2 f q‘ 


arctan 


— arctan 


2 / q 


= 2cro ^ arcsin g' - 
and (2) for g > 1, we have 


= u|(d) 


^( 1 -bg ^cot^il) 


■ COS1? 


In^ 


2 q + q 


q-q' 


-In^ 


q_ 
q' 

2 / (? + g' COS d 


q — q' COS d 


(7) 

( 8 ) 


(9) 


= al^\n 

q 


q + q' 




( 10 ) 


As above, throughout we use the notation g' = |1 — 

We calculate the dynamical normalisations for the 
oblate and the prolate cases using the above equations (7), 
(8), (9) and (10). 

a. The oblate case (q < 1 )—The radial component cr^ and 
the azimuthal component a'^ of the second velocity moments 
in the cylindrical coordinates can be used to obtain a Carte¬ 
sian component that is perpendicular to the symmetry 
axis, namely we have 

2 22 ./ 2.2 /ii\ 

a'x=^RCos (p + sin p. (11) 


The mass-weighted average of a second velocity moment is 
defined by 


r , 2 - ^<1^ \ f ( 12 ) 

I p(r-‘ sin vardr/dp) J 

Using the above definition (equation 12), we obtain for the 
mass-weighted averages of the X and the Z components 
respectively in units of ctq 


wr{q) 

and 

H^z"'(g) 


2/2 H • / 

= — arcsin q 


1 g / g 

=-arctan — 

2 g' 


dy 


sinil 


X ( arctan^ arctan^ 


(13) 


(14) 


Given the mass-weighted averages of the two Cartesian com¬ 
ponents, the line-of-sight component of the second velocity 
moments when the galaxy is viewed at an inclination angle 
of i, with i = 7r/2 defining the edge-on view, is given by 

u-fos)*) = uo[VFz’"'(g) cos^ i + W$^{q) sin^ i], (15) 


For a given axial ratio / of the projected surface density 
given by equation (1), there are infinitely many combina¬ 
tions of inclination angle i and intrinsic axial ratio g that 
are consistent with the projected axial ratio /. Averaging 
equation (15) over this ensemble of models, we have 

f-IT/2 

(u-L) = (7r/2 - imin)“W crL(*)d*, (16) 

•fimin 

where imin is the minimum possible inclination angle for the 
given projected axial ratio /. From the relation 

j = cos ^ -I- g sm I = 1 — g sm ^ (17) 

(for g < 1), we have 

imin = arcsin/', (18) 

where and throughout we use the notation /' = (1 — 
j 2 )i/ 2 . jjitroducing the variable t = sini, we can write the 
inclination-averaged line-of-sight velocity second moment as 

(c^L))/') = <To(7r/2 - arcsin/')"^ 


X / {[wr{q{t))-wr{q{t))r 


+Wf\qm 


dt 




(19) 


The projected surface density for the oblate case of the ax¬ 
isymmetric mass distribution given by equation (6) can be 
written as 


2 

v^obl ^0 


1 


2G ^+ Py"^ 


= ^(g) 


2G ^Py'^ 


-( 20 ) 


where (g) is the inclination averaged value of the axial ratio, 
which can be shown to be 

{q) = (<?)(/) = (7r/2-arcsin/')"^^/V(^l,i;2;/^^ , (21) 

where F is a hypergeometric function. Finally, compar¬ 
ing equation (20) with equation (1) in conjunction with 
the relations given by equations (19) and (21), we obtain 
for the dynamical normalisation of the singular isothermal 
oblate spheroid with a distribution function of the form 
T = T{E,L.) 


A°"'(/) = 


TT 3/2 




{[Wf\q{t))-Wt\qmt^ 


Uf 


+wt\qm 


VT^ 


( 22 ) 


b. The prolate case (q > 1 )— A procedure similar to the 
above procedure for the oblate case can be followed using 
equations (9) and (10) to obtain the dynamical normalisa¬ 
tion for the prolate spheroid. We point out the differences 
in the procedures between the oblate and the prolate cases 
and then list the calculated results. First, the projected ax¬ 
ial ratio /(< 1) is related to the intrinsic axial ratio g(> 1) 


pi / 2 • I 2*2 -N—1 /I I /2 • 2 -N—1 

/ = (cos 1 + q sin i) = (1 -|- g sm i) 


(23) 


Second, the projected surface density for the prolate case of 
the axisymmetric mass distribution given by equation (6) 
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can be written using equation (23) as 
1 


= si fq 


‘2G 


Finally, the requirement of the positivity of the second ve¬ 
locity moments shows that the intrinsic axial ratio q cannot 
be greater than gmax = 3.46717... and thus the minimum 
inclination angle imin for the given projected axial ratio / is 
given by 


= arcsin tmin = arcsin 


/' 


\/ 9max 1 -b 


(25) 


Unfortunately, even for q < Qmax, non-negativity of the dis¬ 
tribution function is not guaranteed. In other words, our 
adopted (simple-minded) dynamical model for the singu¬ 
lar isothermal prolate spheroid has an unphysical feature. 
Nonetheless, our adopted dynamical model was sometimes 
used in the literature as a simple approximation. Further¬ 
more, the construction of a physical dynamical model in con¬ 
junction with a prolate mass model that would mimic the 
SIE is beyond the scope of this paper. Bearing the above in 
mind, we now give the calculated results. The dynamical 
normalisation of the singular isothermal prolate spheroid 
with a distribution function of the form T = J-{E,LP) is 
given by 




-min 

cl 


^2 J2 J/2 


/(1-i^) t 

{[wr\q{t)) - 


+wr\qm 


dt 




where 


W 


proi/ \_19 q q 


‘(g) = ±:iin: 

2q' q-q' 


■ (26) 


(27) 


and 

wr\q) = 


4(7' V q-q’ 


dd- 


sin 


In 


q + q' 


- In 


q + q' cos g 
q — q' cos g 


(28) 


Fig. 1 shows the functional behaviour of the dynamical 
normalisation squared for the following three cases; A(/) = 
A°‘’'(/), A(/) = 0.5A°‘''(/)+0.5AP-'(/), and A(/) = Ap-'(/), 
where A°''*(/) and A‘'^°*(/) are respectively given by equa¬ 
tions (22) and (26) (see section 2.1.1). 


2.1.2 Galaxy luminosity functions, magnification biases, 
and lensing probabilities 

We assume that galaxies of a given morphological type are 
uniformly distributed in comoving space and distributed in 
luminosity following a Schechter (1976) luminosity function 
and that the luminosity of a galaxy has a power-law relation 



Figure 1. Functional behaviour of the lens ‘dynamical normali¬ 
sation’ squared. We consider three cases of the dynamical normal¬ 
ization, namely A(/) = A°'''(/), A(/) = 0.5A°'''(/) -|-0.5AP‘°'(/), 
and A(/) = AP'’°'(/), where A°'''(/) and AP’’°'(/) are respectively 
the dynamical normalisations for the oblate and the prolate scale- 
free isothermal spheroids. 


to its line-of-sight velocity dispersion (i.e. the Faber-Jackson 
relation for early-type galaxies and the Tully-Fisher rela¬ 
tion for late-type galaxies). The differential probability for 
a source with redshift Zs and flux density to be multiply- 
imaged with image multiplicity m and image separation A9 
to A9 + d{A9), due to a distribution of intervening galaxies 
of morphological type g (where g ~ e for early-type galaxies 
and s for late-type galaxies) at redshifts z to z + dz modelled 
by SIEs and described by a Schechter luminosity function, 
is given by 


d^pf}{z, A6>; Zs, SA 
dzd(A9) 


Sm(/)87r^7i?i 


xn*(2:)(l + Z) 


aSz) 


dl 

\b{Q,z)b(z,Zs)'\ 

dz 

o' 

_1 


1 


A9Sz) 
X exp 


A9 


A9Sz) 

A9 

A9Sz] 


1 (Qi7-|-7-|-2)/2 


7/2' 


^ Em (.Zs, Si ,) , 


(29) 


where Sm(/) (hereafter, m = 2,4) is the cross section in 
units of rf (equation 2) for multiple imaging with image 
multiplicity m (section 2.1.1), A9*{z) is the characteristic 
image separation (see below), and Bm(zs, Sv) is the mag¬ 
nification bias (see below). In equation (29), \di/dz\, where 
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£ = //qX( proper time), is given by 
^ = (l + 2 )-'[ 0 ,,(l + 2 )® 

+ (1 — fim — nA)(l + 2 ) + f^A] ^ (30) 


for a cosmology with a classical cosmological constant. In 
equation (29), we have also used the following notations: 
the Schechter (1976) function 


dn{L, z) 

IL 


nt{z) 

U{z) 


L 

U{z) 


exp[-L/L*(z)l, 


(31) 


where a is the faint-end slope, Lt{z) is the characteristic 
luminosity, and nt,{z) is the characteristic comoving number 
density which can be written as 


n*(«) = nt^oe„(z). 


(32) 


The luminosity of a galaxy is assumed to be correlated with 
its line-of-sight velocity dispersion by 


L 

L4z) 


^q0.4[M,(z)-M] 



(33) 


where M^{z) and a*(a) are respectively the characteristic 
absolute magnitude and the characteristic velocity disper¬ 
sion corresponding to the characteristic luminosity L*{z), 
and 7 is the Faber-Jackson or Tully-Fisher exponent. We 
write the characteristic velocity dispersion as 


at(z) = atfiCviz). 


(34) 


In equation (29), the characteristic image separation A9*{z) 
by an L* galaxy is given by 


Ae^z) = 

D{0,zs) 



(35) 


where A(/) is the dynamical normalisation. Notice that the 
characteristic image separation A6t{z) (equation 35) scales 
linearly with the dynamical normalisation A(/), while the 
multiple imaging cross sections Sm{f) scale quadratically 
with A(/) (section 2.1.1). In equations (32) and (34), param¬ 
eters en{z) and e„( 2 ) respectively represent the evolutions of 
the number density and the velocity dispersion of galaxies. 

The magnification bias Bm{zs,Sv) in equation (29), 
namely the factor by which the multiply-imaged sources 
are overrepresented compared with the unlensed sources in 
a flux limited sample because the multiply-imaged sources 
come from intrinsically fainter populations compared with 
the unlensed ones of the same apparent brightness (see, e.g., 
Maoz & Rix 1993; King & Browne 1996), is given by 


section with image multiplicity m that have total magni¬ 
fications greater than /im, for a given lens model. Hence 
\dPm{> ^J,m)/djirnl is the magnification probability distribu¬ 
tion above a minimum possible magnification and it 

follows that Pmi> Pm, min ) = 1. For example, for the SIS 
the magnification probability distribution for double imag¬ 
ing with any magnification ratio of the two images is given 
by |dP®^®(> fJ, 2 )/dfi 2 \ = 8 /p 2 with the minimum total mag¬ 
nification of p 2 ,min = 2. The minimum and maximum total 
magnifications pm,min and pm,max in equation (36) depend 
on observation characteristics as well as the lens model. For 
multiple imaging of a point source, pmax ^ cxd. For the 
CLASS statistical sample, the observational lower limit on 
the ratio of the flux densities of the fainter to the brighter 
image for the doubly-imaged systems is TZmin ~ 0.1. Given 
such an observational limit TT-min, the minimum total mag¬ 
nification for double imaging for the SIS model is given by 

pi:in = 2i±|^. (37) 

-L /Vmin 


For the SIE (i.e. / 7 ^ 1 in equation 1), we numerically 
calculate magnification probability distributions and mini¬ 
mum total magnifications. We do this by solving the sim¬ 
plified lens equation for the SIE (Kormann et al. 1994) 
for millions of source positions randomly distributed on a 
source region that includes the region enclosed by the caus¬ 
tics for each value of the axis ratio /. Fig. 2(a) shows 
several examples of the magnification probability distribu¬ 
tion |dPm(> pm)/dpm| for double and quadruple imaging. 
Fig. 2(b) shows the theoretical minimum total magnifica¬ 
tions for double and quadruple imaging, i.e. for any magnifi¬ 
cation ratios, and the minimum total magnifications for dou¬ 
ble imaging with the CLASS observational limit TZ > 0.1. 

The differential probability for a source to be multiply- 
imaged with image multiplicity m and image separation A9 
to A9 -\- d{A9) due to all the galaxies of type g from z = 0 
to Za, can be obtained from equation (29): 


dpi^\A9-Zs,S4 _ d^p^4}{z,A9\Za,Si,) j 

d{A9) ^ dzd(A9) ^ ’ 

Likewise, the differential probability for a source to be 
multiply-imaged with image multiplicity m and with any 
image separation > A^min® due to galaxies of type g at red- 
shifts z to z + dz is given by 


dpijl\z- Zs, Sy) 

dz 


Sm{f)l67V^r 




R 


3 

H 


xnt{z){l + z)^ 


cra{z) 

c 


Pmi^Zs^SiX) — 


/“Mm,max 

dNzsi^^ Si/jf-lm'] 

1 dPrni^^ /J-m) 


dl 

P(0, z)D{z, Zs) 


dS^ 

djlm 

X 

|_ ^m,min 

1 \ ^ X 1 

-1 

dz 

D{0,Za) 


X- dpLrt 


dN.4> S4 


dS, 


(36) 


X Bm (^Zs , iS'fy ) 

/•AVin d?p4\z,A9-Za,S4 


(m = 2, 4). Here Nza{> S4 is the intrinsic number- 
flux density relation, i.e. the integrated source counts as 
a function of flux density S^, for the source population 
at redshift Zs- Parameter prn denotes the total magnifica¬ 
tion for multiple imaging with image multiplicity m, namely 
the sum of the absolute magnifications of the m images. 
Pm(> pm) is the fraction of the multiple imaging cross 
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dzd{A9) 


d{A9), 

(39) 


^ An upper limit on A0 is not necessary since the lensing prob¬ 
ability due to galaxy lenses (given by equation 29) is negligibly 
small beyond an observational limit, e.g. several arcseconds. 
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(a) Magnification probability distributions 



(b) Minimum total magnifications 



Figure 2. (a) Magnification probability distributions (MPDs). The six curves from right to left respectively correspond to MPDs for 
quadruple imaging by e = 0.05, 0.3 and 0.6 SIEs and for double imaging by e = 0.05, 0.3 and 0.6 SIEs. The thin solid lines representing 
/^2 behaviours and the thin dashed lines representing behaviours fit respectively well the numerical data points for double 

imaging and quadruple imaging at high magnifications in agreement with the expectations from analytical calculations, (b) Minimum 
total magnifications for double and quadruple imaging by SIEs as a function of ellipticity e. The solid line shows the minimum total 
magnifications for double imaging for any magnification ratios between the two images. The dotted line shows the minimum total 
magnifications for double imaging only for the fainter image to the brighter image magnification ratio > 0.1, which is the observational 
limit for the final CLASS statistical sample. The dashed line shows the minimum total magnifications for quadruple imaging for any 
magnification ratios between the four images. 


where the integrand in the last piece of the right-hand side 
is given by equation (29). Finally, using equation (39) we 
obtain the integrated probability for the multiple imaging 
with image multiplicity m due to galaxies of type g\ 


S9)(^ Q 

-'m ? ‘Ji/ 


dpif’{z-, Zb, S„ 

dz 


-dz. 


(40) 


2.1.3 Summary of model parameters 

The lensing probabilities given by equations (29), (38), (39) 
and (40) are for a population of galaxies of a specific mor¬ 
phological type that are assumed to be described by a set of 
Schechter parameters n*, a, and M* with corresponding L* 
and (7* along with the Faber-Jackson or Tully-Fischer expo¬ 
nent 7 (equation 33). We assume two populations of galax¬ 
ies separately contributing to multiple imaging, namely the 


early-type (i.e. ellipticals and SOs) and the late-type (i.e. 
spirals and other non-E/SOs) populations. Separately tak¬ 
ing into account the two populations is required. This is 
because the early-type and the late-type galaxies are dy¬ 
namically different, satisfy respectively the Faber-Jackson 
and the Tully-Fisher relations, and have different character¬ 
istics in multiple imaging: the multiple imaging rate by the 
early-type population is « 3-6 times higher and the early- 
type galaxies produce on average larger image splittings. 


The parameters of our statistical lensing model can be 
summarized as follows: (1) up to three cosmological param¬ 
eters: • matter density flm, • density of a classical cosmo¬ 
logical constant Ua or dark energy density fla, with • its 
constant equation of state w = Pa, (pressure)/pa, (energy den- 
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sity);® (2) six parameters for the early-type galaxy popu¬ 
lation, i.e. early-type Schechter parameters and 

Faber-Jackson exponent ypj, early-type mean mass 
ellipticity and relative frequency of the oblatenesses 
among early-type galaxies; (3) six parameters for the 
late-type galaxy population, i.e. late-type Schechter param¬ 
eters and Tully-Fisher exponent 7tf, late- 

type mean mass ellipticity , and relative frequency of the 
oblatenesses V^i among late-type galaxy halos.^ 

In addition to the above parameters, our statistical leas¬ 
ing model includes the parameters that can be used to model 
certain aspects of galaxy evolutions, i.e., and ei^\z) 

where g = e,s (see equations 32 and 34). In this work, we 
assume no evolution of the comoving number densities of 
galaxies for the redshift range which is relevant to lensing 
by CLASS sources, i.e., en\z) = 1 for 0 < 2<1. Here 
note that any significant evolution of the comoving num¬ 
ber density of late-type galaxies has little effect on the de¬ 
rived values of cosmological parameters. Note also that while 
a significant evolution of the comoving number density of 
early-type galaxies would have a significant effect on the 
derived values of cosmological parameters, (our interpreta¬ 
tion of) current observational evidence suggests that there 
has been little evolution of the comoving number density 
of early-type galaxies since « ~ 1 (see the last paragraph 
of section 3.3). However, evolution of the comoving number 
density of early-type galaxies is a main source of uncertainty 
and its systematic effects on the derived values of cosmo¬ 
logical parameters are discussed in section 5.3.3. Since the 
characteristic velocity dispersions will be treated as free pa¬ 
rameters (section 4) and lens redshifts are intermediate, it 
should be understood that parameters ei^\z) are irrelevant 
and we are constraining parameters for 0.3<2< 1. Fi¬ 
nally, in addition to the cosmological parameters and the 
lens parameters, there are parameters used to describe the 
CLASS source population, namely the redshift distribution 
and the number-flux-density relation, which can be found 
in section 3.2. 


2.2 Method of statistical analysis 

For a statistical sample that contains Al multiply-imaged 
sources and Nu unlensed sources, the likelihood £ of the ob¬ 
servation data given the statistical lensing model is defined 


by 

Nu 

P 

-i 

ln£ 

- El- 

k=i 

1 - 

pi^“\zs,k,Sn,k) 

m=2,4 


® Here all the cosmological densities refer to the fractions of the 
present critical density. 

^ In this work, we are considering a single-component galaxy 
model. This means that for a spiral galaxy the mass model is 
intended to approximate the total mass distribution by the halo, 
the bulge and the disk of the galaxy. Consequently, the model 
parameters (i.e., velocity dispersion and projected ellipticity) are 
effective quantities corresponding to the combined lensing effects 
by the components of the galaxy. Discussions on the lensing ef¬ 
fects by the components of a spiral galaxy can be found, e.g., in 
Keeton & Kochanek (1998). 


JVl 

+ ^ In Aer, Zs,i. S,,il (41) 

l = \ 

where X/m =2 S^^k) fs the sum of the double and 

quadruple imaging probabilities (equation 40) for the fc-th 
unlensed source, and Spm’^’^\{zi), A9i; Zs,i, S,y^i] is the dif¬ 
ferential probability (equation 29, 38, or 39) for the Lth 
multiply-imaged source, which is a specific probability in¬ 
cluding the observed configuration of the source, i.e. its im¬ 
age separation AOi, its image multiplicity mi (and its lens 
redshift zi if it is known). 

In equation (41), for the unlensed sources the proba¬ 
bility to be used is the sum of the probabilities due to the 
early-type and the late-type populations, namely 

pin ^{Zs,k^Sn,k) = pin {Zs ^k ^ Sn,k) pin {Zs ^k ^ Si/^k\ 

= 1,2, ...,Au. (42) 

For the multiply-imaged sources, the suitable differential 
probability can be given by 

[{Zl), A9i; Zs,l, Sn,l] 

Spi^l [{Zl), A9l-Zs,l, Sn,l], 

/ = 1,2,...,Al, (43) 

where parameters and are the weights given to 
the early-type and the late-type populations, respectively, 
satisfying = 1. If the lensing galaxy type is known 

to be an early-type (late-type), = 1 = 1). If the 

lensing galaxy type is unknown, we use + 

^pini\ {g = e, s), if the lens redshift is known, whereas we 
use = 0.8 and = 0.2 otherwise. 

From equation (41), a “x^” is defined by 

X^ = —2ln£ [-l-Gaussian prior on (44) 

As indicated in the bracket in equation (44), the x^ "'iU 
sometimes include a Gaussian prior probability term of 
the late-type characteristic velocity dispersion ai'’ . We de¬ 
termine ‘best-fit’ model parameters by minimizing the x^ 
(equation 44) and obtain confidence limits on the model pa¬ 
rameters using the usual Ax^ (= X^ ~ Xmin)“Static, where 
Xmin is the global minimum x^ for the best-fit parameters. 


3 INPUT DATA 

In this section, we summarize data both from CLASS obser¬ 
vations and recent large-scale observations of galaxies that 
will be used as the input for the statistical lensing model 
described in section 2. In section 3.1, we briefly review the 
definition of the final CLASS statistical sample that satisfy 
well-defined observational selection criteria (Browne et al. 
2003) and summarize the properties of the multiply-imaged 
sources in the sample. In section 3.2, we summarize observa¬ 
tional information on flat-spectrum radio sources including 
that from CLASS observations and estimate the distribution 
of CLASS radio sources in redshift and flux density. In sec¬ 
tion 3.3, we survey recent observational results for galaxy 
luminosity functions and obtain estimates of galaxy lumi- 
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nosity functions per morphological type based on the data 
from the literature. 


3.1 The final CLASS statistical sample 

The Cosmic Lens All-Sky Survey is the largest completed 
(radio-selected) galactic mass scale gravitational lens search 
project. The CLASS project along with its predecessor 
project the Jodrell-Bank-VLA Astrometric Survey (JVAS; 
see, e.g., King et al. 1999) identified 22 multiply-imaged sys¬ 
tems out of a total of 16521 radio sources targeted (Myers 
et al. 2003; Browne et al. 2003). Out of the entire CLASS 
sample including the JVAS sources, a subsample of 8958 
sources containing 13 multiply-imaged systems constitutes 
a statistical sample that satisfies well-defined observational 
selection criteria (Browne et al. 2003). The final CLASS sta¬ 
tistical sample is the largest sample that can be used for the 
purpose of statistical analyses of gravitational lensing. The 
reader is referred to Myers et al. (2003) and Browne et al. 
(2003) for the description of the CLASS observations includ¬ 
ing the selection of the targets, the gravitational lens can¬ 
didate selection process, the follow-up observations of the 
lens candidate sources, and the properties of the finally con¬ 
firmed individual multiply-imaged systems. From Browne 
et al. (2003), the final CLASS statistical sample satisfies the 
following well-defined selection criteria that are observation- 
ally reliable: 

(i) The spectral index between 1.4 GHz and 5 GHz is 
flatter than —0.5, i.e. a > —0.5 with oc u°‘. 

(ii) The total flux density of (the components of) each 
source > 30 mjy at 5 GHz. 

(iii) The total flux density of (the components of) each 
source > 20 mJy at 8.4 GHz. 

(iv) For multiply-imaged systems, (a) the compact radio¬ 
core images have separations > 300 milli-arcseconds, and 
(b) for doubly-imaged systems the ratio of the flux densities 
of the fainter to the brighter images is > 0.1. 

The final CLASS statistical sample includes JVAS sources 
which essentially form the bright tail with their flux densities 
> 200 mJy at 5 GHz. However, not all JVAS sources are 
included in the final CLASS statistical sample because of 
the strict application of the above criteria (i) and (iv). 

The properties of the 13 multiply-imaged systems con¬ 
tained in the final CLASS statistical sample are summa¬ 
rized in Table 1. Notice that two out of the four multiply- 
imaged sources in the JVAS (statistical) sample used by 
Helbig et al. (1999), namely 0414-1-054 and 1030-1-074, are 
excluded because 0414-1-054 has too steep a spectral in¬ 
dex and 1030-1-074 has two images whose fainter-to-brighter 
image flux-density ratio is less than 0.1. Notice also that 
JVAS system 2114-1-022 which was excluded by Helbig et 
al. (1999) is included along with other multiple-galaxy lens 
systems. The total flux densities, the lens redshifts (if mea¬ 
sured), the source redshifts (see below for sources without 
measured redshifts), the image separations, the image mul¬ 
tiplicities, and the lensing galaxy types (if determined), as 
given in Table 1, are all used (through the likelihood func¬ 
tion defined in section 2.2) to constrain the parameters of 
the statistical lensing model. However, the measured image 
separations of 1359-1-154, 1608-1-656, and 2114-1-022 are not 


used because the observed angular sizes are due to multi¬ 
ple galaxies within their critical radii and thus should not 
be used to constrain parameters pertaining to single galax¬ 
ies. Notice particularly that lens modelling of these systems 
shows that the inferred critical radii of the secondary (and 
tertiary) galaxies are comparable in size to those of the cor¬ 
responding primary galaxies; the ratio of the smallest to 
the largest critical radii ranges from 66% to 81% assum¬ 
ing singular isothermal mass models for individual galax¬ 
ies [see Rusin et al. (2001a) for 1349-1-154, Koopmans & 
Fassnacht (1999) for 1608-1-656, and Chae, Mao, & Augusto 
(2001) for 2114-1-022]. Furthermore, the inferred critical radii 
of the primary galaxies for the multiple-galaxy lens sys¬ 
tems are similar to those of the early-type lensing galaxies 
for the single-galaxy lens systems for which no compara¬ 
bly massive secondary galaxies are observed or required for 
lens modelling; specifically, the average critical radius of the 
primary galaxies for 1349-1-154, 1608-1-656, and 2114-1-022 
is 1.34 arcsec while that for eight (presumably) early-type 
single-galaxy lens systems, namely 0445-1-123, 0631-1-519, 
0712-b472, 0850-b054, 1152-bl99, 1422-b231, 1933-b503, and 
2319+051 (see below why 2045+265 is excluded), is 1.23 
arcsec. These results suggest that discarding the observed 
image separations for the multiple-lens systems 1349+154, 
1608+656, and 2114+022 would be an appropriate way of 
interpreting the data. Likewise, the image multiplicities of 
the above three systems are not used to constrain mean 
mass ellipticity of galaxies because asymmetries in the to¬ 
tal potentials of these systems are caused by combinations 
of multiple potentials. The image separation of 2045+265 is 
not used because of present uncertainties in the details of 
the lensing scenario. For the multiply-imaged sources with¬ 
out measured source redshifts, we take Zs = 2 which is the 
mean source redshift for the multiply-imaged sources with 
measured source redshifts. 

3.2 Global properties of flat-spectrum radio 
sources 

As seen in section 2.1.2, the calculation of the lensing prob¬ 
ability (equations 29, 38, 39, or 40) for a source requires 
the knowledge of the differential number-flux-density rela¬ 
tion {\dNz^{> Sv)/dS,z\) that the source follows and the 
redshift of the source (zs). More specifically, the differential 
number-flux-density relation is required in the calculation 
of the magnification bias (equation 36). In particular, the 
differential number-flux-density relation needs to be known 
from the flux density of the source to lower flux densities 
(formally to the zero flux density). Furthermore, the dif¬ 
ferential number-flux-density relation needs to be known 
as a function of source redshift. However, systematic mea¬ 
surements have been done only for a few dozens of CLASS 
sources (Marlow et al. 2000). For this reason, we are led to 
ignore any redshift dependence of the differential number- 
flux-density relation. We also ignore weak dependences of 
the differential number-flux-density relation on the spectral 
index of the source for the sake of simplicity. 

The final CLASS statistical sample is well-described 
at u = 5 GHz by |dA(> S 5 )/dS 5 \ oc (S’s/S'g)"’’ with 
rj = 2.07 ± 0.02 for S'® < S's < 1 Jy where S'® = 30 mJy 
(McKean et ah, in preparation). For flux densities lower 
than Ss, McKean et al. (in preparation) find rj = 1.97+0.14 
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Table 1. Multiply-imaged systems in the final CLASS statistical sample. The ratio of the flux densities of the images is given for the 
doubles only. Probable lens galaxy morphological-type identification is given using the codes, ‘e’ for an early-type and ‘s’ for a spiral. 
The references for the morphological type identifications are: (1) Browne et al. (1993); (2) Fassnacht & Cohen (1998); (3) Impey et al. 
(1996); (4) Myers et al. (1995); (5) Sykes et al. (1998); (6) Augusto et al. (2001); (7) Rusin et al. (2001b). 


Source 

Survey 

Total flux 
density (mJy) 

Source 
redshift Zs 

Lens 

redshift zi 

Maximum image 
separation (arcsec) 

Image flux- 
density ratio 

Number of 
images 

Lensing 
galaxy type 

0218-1-357 

JVAS 

1480. 

0.96 

0.68 

0.334 

0.26 

2 

si 

0445-1-123 

CLASS 

50. 

- 

- 

1.33 

0.14 

2 

- 

0631-1-519 

CLASS 

88. 

- 

- 

1.16 

0.15 

2 

- 

0712-1-472 

CLASS 

30. 

1.34 

0.41 

1.27 

- 

4 

e2 

0850-1-054 

CLASS 

68. 

- 

- 

0.68 

0.14 

2 

- 

1152-1-199 

CLASS 

76. 

1.019 

0.439 

1.56 

0.33 

2 

- 

1359-1-154 

CLASS 

66. 

3.235 

- 

1.65 

- 

6 

- (3 Gs) 

1422-1-231 

JVAS 

548. 

3.62 

0.34 

1.28 

- 

4 

e3 

1608-1-656 

CLASS 

88. 

1.39 

0.64 

2.08 

- 

4 

e (2 Gs)4 

1933-1-503 

CLASS 

63. 

2.62 

0.755 

1.17 

- 

4 

e®? 

2045-1-265 

CLASS 

55. 

- 

0.867 

1.86 

- 

4 

- 

2114-1-022 

JVAS 

224. 

- 

0.32/0.59 

2.57 

0.33 

2? 

e (2 Gs)® 

2319-1-051 

CLASS 

76. 

- 

0.624/0.588 

1.36 

0.18 

2 



based on a Very Large Array (VLA) mapping of selected 
regions of the sky down to = 5 mjy. The differential 
number-flux-density relation can also be estimated from 
model flat-spectrum radio luminosity functions found in the 
literature. Dunlop & Peacock (1990, hereafter DP90) derived 
flat-spectrum radio luminosity functions based on redshift 
measurements of hundreds of flat spectrum radio sources 
with S 2.7 > 100 mJy. Specifically, DP90 presented five free¬ 
form radio luminosity functions and two parametric radio 
luminosity functions (i.e. a pure luminosity-evolution model 
and a model with both luminosity evolution and density evo¬ 
lution). Waddington et al. (2001) tested the DP90 models 
against the redshift data for a Leiden Berkeley Deep Sur¬ 
vey (LBDS) Hercules sample with Si ,4 > 1 mJy. Wadding¬ 
ton et al. ( 2001 ) find that the measured redshifts for the 
LBDS Hercules sample are consistent with two of the DP90 
models, namely free-form model number 4 (FF-4) and free¬ 
form model number 5 (FF-5). Fig. 3 displays the measure¬ 
ments by McKean et al. (in preparation) and the predic¬ 
tions by the DP90 FF-4 and FF-5 models for the differential 
number-fiux-density relation for flat-spectrum radio sources 
for S 5 > 1 mJy. As shown in Fig. 3, for flux densities below 
30 mJy the measurements by McKean et al. (in preparation) 
are roughly consistent with the extrapolation by the DP90 
FF-5 model but inconsistent with the DP90 FF-4 model. 

Independent pieces of available evidence that can be 
used to estimate the redshift distribution of the CLASS 
sources are concordant. Based on spectroscopic observations 
of a representative subsample of 42 CLASS sources with 
flux densities from Ss = 25-50 mJy, Marlow et al. (2000) 
obtained a mean redshift (z) — 1.27 with a dispersion of 
0.95, at a completeness level of 64%. Since the DP90 FF- 
4 and FF-5 models are consistent with the LBDS redshift 
data with > 1 mJy (see above), we are allowed to use the 
DP90 FF-4 and FF-5 models to infer mean redshifts for flat- 
spectrum sources with Ss ^ 1 mJy. Fig. 4 shows the mean 
redshift for fiat-spectrum sources as a function of flux den¬ 
sity predicted by the DP90 FF-4 and FF-5 models. Fig. 4 
also shows results from various existing redshift observations 
of flat-spectrum radio sources, including the Marlow et al. 



Figure 3. Differential number—flux-density relation \dN{> 
S:i)/dSs\ for fiat-spectrum radio sources. The measurement by 
McKean et al. (in preparation) is compared with the predictions 
by DP90’s FF-4 and FF-5 models (see section 3.2). 


(2000) result and the results by Falco, Kochanek, & Munoz 
(1998), Henstock et al. (1997), and Jackson et al. (2002). 
From Fig. 4 we see that the prediction by the DP90 FF-5 
model agrees well with the independent data while the DP90 
FF-4 model is marginally inconsistent with the data, simi- 
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Figure 4. Mean redshifts of flat-spectrum radio sources as a 
function of flux density. The measured mean and its Icr error 
from various observations are shown using a fllled dot (CLASS 
subsample; Marlow et al. 2000), fllled triangles (Falco, Kochanek, 
& Munoz 1998), filled squares (Parkes quarter-Jansky sample; 
Jackson et al. 2002), and an open square (the second Caltech- 
Jodrell Bank sample; Henstock et al. 1997). The predictions by 
DP90’s FF-4 and FF-5 models (see section 3.2) are respectively 
shown using dotted and solid lines. 


larly to the case for the differential number-flux-density re¬ 
lation (see above). This agreement increases our confidence 
in the Marlow et al. (2000) measured mean redshift for the 
CLASS sample. Therefore, we adopt the measurement by 
Marlow et al. (2000) for a mean redshift for the CLASS 
sources; i.e., we take a mean 2 = 1.27 ±0.18 where the esti¬ 
mated error is from the measured dispersion 0.95. 

Given that there are only 26 unbiased CLASS sources 
that have measured spectroscopic redshifts (Marlow et al. 
2000), the distribution of the redshifts of the CLASS sources 
is poorly known. However, since mean redshifts for flat- 
spectrum sources are nearly the same for various samples 
regardless of their flux density ranges, as shown in Fig. 4, 
we may combine the samples assuming that the redshift dis¬ 
tributions of the samples are similar. The combined sam¬ 
ple comprising the samples shown in Fig. 4 has 747 flat- 
spectrum source redshifts with S 5 > 30 mjy. A histogram 
of the 747 redshifts is shown in Fig. 5. Figure 5 also shows 
redshift distributions predicted by the DP90 FF-4 and FF- 
5 models as well as a Gaussian model whose peak redshift 
and width are determined from the mean and dispersion 
measured by Marlow et al. (2000). From Fig. 5 we see that 
the Gaussian model adequately describes the redshift dis¬ 
tribution for the combined sample. Therefore, we adopt the 
Gaussian model to approximate the redshift distribution of 
the GLASS sources. 



Figure 5. Redshift distribution of flat-spectrum sources. Red¬ 
shift distributions predicted by the DP90 FF-4 and FF-5 models 
given for S 5 = 50 mJy and a Gaussian model distribution (see 
section 3.2) are overplotted on a histogram of 747 flat-spectrum 
source redshifts from various observations (Henstock et al. 1997; 
Falco, Kochanek, & Muiioz 1998; Marlow et al. 2000; Jackson et 
al. 2002). 


3.3 Galaxy luminosity functions per 

morphological type and their evolution since 

2 ~ 1 

Use of reliable galaxy luminosity functions (LFs) is vital in 
analyses of statistical lensing. For example, use of an un¬ 
derestimated (overestimated) comoving number density of 
galaxies would mislead us to overestimate (underestimate) a 
dark energy density. In particular, given that galaxy popula¬ 
tions of different morphological types contribute to multiple 
imaging in distinctively different ways (see section 2.1.3), 
it is required that type-specific galaxy LFs are reliably de¬ 
termined. For this reason, in this subsection we thoroughly 
review recent observational results on galaxy LFs and ob¬ 
tain estimates of the early-type and the late-type LFs based 
on the data from the recent literature. We consider first the 
estimation of the local type-specific LFs and then review 
recent observational studies on the evolution of early-type 
galaxies since 2 ~ 1. 

We consider two approaches for estimating the local 
type-speciflc LFs. The first approach is to estimate the lo¬ 
cal total LF and decompose it into an early-type LF and a 
late-type LF. The second approach is to simply adopt type- 
specific LFs measured in large-scale galaxy surveys. The ad¬ 
vantage of the first approach is that observational informa¬ 
tion on the total LFs is abundant and now appears to pro¬ 
duce a converging result on the total LF, owing to more 
recent large galaxy redshift surveys particularly including 
the Sloan Digital Sky Survey (SDSS) and the Two Degree 
Field Galaxy Redshift Survey (2dFGRS). However, infer¬ 
ence of the type-specific LFs from the total LF relies on an 
independent knowledge of relative type-specific LFs. By rel¬ 
ative type-specific LFs, we refer to the relative magnitudes 
of the characteristic number densities of the type-specific 
LFs. Thus, for the determination of the relative type-specific 
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LFs, the sample size of the galaxies used is less relevant but 
the accuracy of the galaxy classification and the photometry 
are important. Direct use of the measured type-specific LFs 
from large-scale galaxy surveys in which galaxies were clas¬ 
sified (the second approach) suffers from a relatively small 
number of independent measurements and uncertainties in 
classifications of large numbers of galaxies in each measure¬ 
ment. 

We now consider estimating the total LF taking into ac¬ 
count recent data from the literature. Observational results 
from various past galaxy surveys had systematic differences 
among the estimates of the local total LF. Not only were the 
measured values of the three Schechter parameters M*,o, o 
and n*,o significantly different among the various observa¬ 
tional results, but also the estimates of the integrated lumi¬ 
nosity density given by 

= ^ + 2) (45) 

differed by up to a factor of two.® The past results indi¬ 
cate that, aside from any measurement errors (e.g. photo¬ 
metric errors, missing or misclassification of galaxies), the 
various determinations of the total LF suffered from large- 
scale structures in the (local) Universe. The effects of large- 
scale structures on the determinations of the total LF ap¬ 
pear to have caused overall differences between measured 
LFs based on wide and shallow/sparse surveys (e.g. Love- 
day et al. 1992; Marzke et al. 1998, hereafter M98) and those 
based on narrow and deep surveys (e.g. Lilly et al. 1995; El¬ 
lis et al. 1996; Geller et al. 1997; Zucca et al. 1997), in the 
sense that normalisations of the LF for shallow samples are 
significantly lower than those for deeper samples. Systematic 
differences in the number counts of galaxies between those 
towards the northern Galactic cap and the southern Galactic 
cap also provide evidence of large-scale structures (Yasuda 
et al. 2001; Norberg et al. 2001). Thus, past determinations 
of the total LF depended, to a varying extent, on the coor¬ 
dinates defining the sky patches observed in the surveys as 
well as the angular areas and depths of the surveys. 

However, several recent galaxy surveys contain large 
numbers of galaxies because of the large areas and/or great 
depths of the surveys. In addition, for a few of the recent 
galaxy surveys, CGD-based images of galaxies were used 
to measure the galaxy LFs. More importantly, the results 
from several reliable observations appear to be concordant. 
We consider the 2dFGRS (Cross et al. 2001, Norberg et 
al. 2001), the SDSS (Blanton et al. 2001, Yasuda et al. 
2001), the Century (Brown et al. 2001), the ESP (ESO Slice 
Project; Zucca et al. 1997), and the LCRS (Las Campanas 
Redshift Survey; Lin et al. 1996) surveys. The 2dPGRS 
(Cross et al. 2001, Norberg et al. 2001), the SDSS (Blan¬ 
ton et al. 2001, Yasuda et al. 2001), and the LCRS (Lin 
et al. 1996) observations contain large numbers of galax¬ 
ies (over 10,000 galaxies each). The SDSS (Blanton et al. 
2001; Yasuda et al. 2001), the Century (Brown et al. 2001), 


® Differences in the measurements of each Schechter parameter 
are in part due to degeneracies among the three Schechter param¬ 
eters in fitting the data. However, the calculated quantity given 
by equation (45) would not be significantly affected by such de¬ 
generacies. 


and the LCRS (Lin et al. 1996) observations use CCD-based 
photometry. Furthermore, the SDSS (Blanton et al. 2001; 
Yasuda et al. 2001) and the Century (Brown et al. 2001) 
observations use multiple passbands. All of the above re¬ 
cent surveys except for the LCRS survey (Lin et al. 1996) 
are based on passbands similar to the Johnson-Morgan B 
passband in wavelength coverage. The measured values of 
the Schechter parameters and the luminosity density (equa¬ 
tion 45) from these recent surveys are summarized in Ta¬ 
ble 2. Notice that the characteristic magnitudes for all the 
above surveys except for the LCRS are quoted in the bj pho¬ 
tometric system using the relevant photometric conversions 
quoted in the table. 

The Norberg et al. (2001) result was obtained by nor¬ 
malising their LF to the combined galaxy number counts 
from the 2dFGRS southern and northern Galactic caps (over 
110,500 galaxies). Likewise, the Yasuda et al. (2001) result 
was obtained by correcting the LF calculated by Blanton et 
al. (2001) to match the number counts of galaxies in SDSS 
(~ 900,000 galaxies). The Norberg et al. (2001) result in¬ 
cluded modelling and correcting for galaxy evolutionary ef¬ 
fects, and the SDSS results were derived using Petrosian 
magnitudes of galaxies (see Blanton et al. 2001), which over¬ 
come dependences of isophotal magnitudes on the surface- 
brightness distribution of galaxies. The Norberg et al. (2001) 
and the Yasuda et al. (2001) results are independent of and 
consistent with each other and probably represent present 
best estimates of the local total LF. The results of the Cen¬ 
tury (Brown et al. 2001) and the ESP (Zucca et al. 1997) 
surveys are also consistent with the Norberg et al. (2001) 
and the Yasuda et al. (2001) results; however, their statisti¬ 
cal errors are larger. The result of the LCRS (Lin et al. 1996), 
which is based on an R passband, appears to be significantly 
different from those of the other surveys. We consider the 
LCRS (Lin et al. 1996) because of the large comoving vol¬ 
ume covered by the survey and the fact that it is one of 
the few CCD-based imaging surveys. While the characteris¬ 
tic magnitude and normalisation of the LE depend on the 
wavelength coverage of the passband used, recent determi¬ 
nations of LEs at multiple passbands based on CCD pho¬ 
tometry did not show any evidence for dependence of the 
faint-end slope on passbands of optical wavelengths (Blan¬ 
ton et al. 2001; Brown et al. 2001). Thus, the significantly 
shallower faint-end slope obtained from the LCRS (Lin et 
al. 1996) compared with those from the other surveys given 
in Table 2 might appear to be a significant obstacle to a con¬ 
cordant total LF. However, using the SDSS data, Blanton 
et al. (2001) illustrate dependences of the chosen isophotal 
limits on the derived LF parameters. In particular, Blan¬ 
ton et al. (2001) argue that the result of Lin et al. (1996) 
was significantly affected by their chosen relatively bright 
isophotal limit. We use the Norberg et al. (2001) result and 
the Yasuda et al. (2001) result to estimate the Schechter pa¬ 
rameters of the local total LF. Although the errors of the 
three Schechter parameters may be correlated, we ignore 
such correlations and calculate weighted means of the val¬ 
ues of the parameters obtained by Norberg et al. (2001) and 
Yasuda et al. (2001). This result is our estimate of the local 
total LF given in Table 2.® 


® We notice that our estimate of the Schechter parameters is also 
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Table 2. Summary of recent measurements of the local total luminosity function (1. Cross et al. 2001; 2. Norberg et al. 2001; 3. Blanton 
et al. 2001; 4. Yasuda et al. 2001; 5. Brown et al. 2001; 6. Zucca et al. 1997; 7. Lin et al. 1996) 


# 

Survey 

A7*,o - SlogiQ h 

a 

(10-2 /i3 Mpc-3) 

jo 

(10^ h Lq Mpc-3) 

1 

2dFGRS 

-19.75 ±0.05 (fej) 

-1.09 ± 0.03 

2.02 ± 0.02 

2.24 ±0.12 

2 

2dFGRS 

-19.66 ±0.07 (bj) 

-1.21 ± 0.02 

1.68 ±0.08 

1.90 ±0.18 

3 

SDSS 

-19.71 ±0.05 (fej)“ 

-1.26 ± 0.05 

2.06 ±0.23 

2.57 ±0.34 

4 

SDSS 

-19.71 ±0.05 (fej)“ 

-1.26 ± 0.05 

1.81 ±0.16 

2.26 ±0.26 

5 

Century 

-19.5 ±0.09 (fej)*’ 

-1.07 ±0.09 

2.0 ±0.3 

1.74 ±0.32 

6 

7 

ESP 

LCRS 

-19.611^8 (bj) 

-20.29 ±0.02 (Rkc) 

r)Q + 0.06 

-0.70 ±0.05 

2.0 ±0.4 

1.9 ±0.1 

2.19 ±0.48 

1.4 ±0.1 


Our estimate 

-19.69 ±0.04 (fej) 

-1.22 ± 0.02 

1.71 ±0.07 

2.02 ±0.15 


“ The SDSS g* magnitude was converted to the bj magnitude using bj = g* + 0.12 + 0.16(B — V) and B — V = 0.94 (Norberg et al. 

2001) 

^ The Century V magnitude was converted to the bj magnitude using bj — R = 1.2 and V — R = 0.53 (Brown et al. 2001). 


To extract type-specific LFs from our estimated total 
LF, we define the ratio of the integrated luminosity density 
of a population of galaxies of type g {g = e, s) to the total 
integrated luminosity density, i.e., 


An) - jo ^ 


"■*,0 j^q0.4(M,.o-M 

n*,o 


+ 2 ) 

T{a + 2) 


(46) 


where the last equality follows from equation (45). Notice 
that the fractional luminosity density given by equation (46) 
is unaffected either by an error in the absolute photometry 
or large-scale galaxy-number-density fluctuations. The accu¬ 
racy in the measurement of the fractional luminosity density 
in a survey is limited mainly by errors in relative photometry 
and errors in the galaxy classification. It follows from equa¬ 
tion (46) that the type-specific normalisation is given 
by 




-M, o) r(a -|- 2) 


r(a(9) -b 2 )’ 


(47) 


Given an estimate of the total LF, independent estimates 
of parameters and along with the the fractional 

luminosity density will allow us to estimate the type- 
specific normalization n$,®Q using equation (47). 

We now consider direct measurements of the early-type 
Schechter parameters from recent galaxy surveys and derive 
corrected early-type Schechter normalisations using equa¬ 
tion (47) from the directly measured early-type Schechter 
parameters and our estimate of the total LF given in Ta¬ 
ble 2. Direct measurements of the early-type LF from galaxy 
surveys are relatively few and have relatively larger uncer¬ 
tainties mainly because of the challenging observational task 
of accurately classifying large numbers of galaxies. The ob¬ 
servational techniques of classifying galaxy types used in 
various galaxy surveys are based on morphological appear¬ 
ances, surface brightness profiles, prototype spectra, spec¬ 
tral energy distributions, and/or colours of galaxies. Iden¬ 
tifying early-type galaxies based on morphological appear¬ 
ances or surface brightness profiles is probably the most re¬ 
liable method, since the morphological appearances and the 


in good agreement with just recently published galaxy counts in 
the Millennium Galaxy Catalogue (Liske et al. 2002). 


fight profiles of spirals, which constitute the most numer¬ 
ous population, are well distinguished from those of early- 
types. We consider the results from the Second Southern Sky 
Redshift Survey (SSRS2; M98) and the 2dFGRS (Folkes et 
al. 1999; Madgwick et al. 2001).^° The SSRS2 (M98) is a 
wide-angle, shallow {z < 0.05) survey of 5404 galaxies in 
the northern and southern Galactic caps that are brighter 
than tub ~ 15.5. The relatively bright magnitude limit used 
for the SSRS2 allowed M98 to classify the galaxies in their 
sample through the visual inspection of morphological ap¬ 
pearances with more than 99% completeness. M98’s deter¬ 
mination of the early-type LF probably represents the best 
estimate of the local early-type LF. However, as they discuss, 
the M98 sample would not provide a fair number density of 
galaxies because the relatively small volume covered by their 
sample suffers from large-scale density fluctuations. In this 
respect, the M98 result would only be valuable as a relative 
early-type LF. M98’s measured values of the local early-type 
LF parameters can be found in Table 3. 

Apart from the M98 study, recent observational studies 
on galaxy LFs per type did not classify galaxies through 
morphological appearances but used classification tech¬ 
niques that are based on photometric and spectroscopic in¬ 
formation on galaxies (e.g. Folkes et al. 1999; Blanton et al. 
2001; Fried et al. 2001; Madgwick et al. 2001; Bernard! et 
al. 2003). We consider the results by Folkes et al. (1999) 
and Madgwick et al. (2001) which are based on the 2dF- 
GRS data. Folkes et al. (1999) and Madgwick et al. (2001) 
have employed the so-called “principal component anal¬ 
ysis (PGA)” of galaxy spectra, which essentially maxi¬ 
mally quantifies differences between the spectra of galax¬ 
ies. Kochanek, Pahre, & Falco (2000) have recently argued 
that the physical properties of apertures (particularly their 


An estimate of the Schechter parameters per galactic type is 
not yet reported from the SDSS, although some classifications 
of galaxies have been considered (Blanton et al. 2001; Bernard! 
et al. 2003). In particular, we note that the sample of early-type 
galaxies used by Bernard! et al. (2003) in their study of early-type 
galaxies were selected too stringently (for example, luminosities 
of galaxies in the sample spread over only 4 magnitudes in the 
brightest range) and thus not suitable to derive a ‘fair’ LF of 
early-type galaxies. 
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Table 3. Summary of estimates of the local early-type luminosity function (1. M98; 2. Folkes et al. 1999; 3. Madgwick et al. 2001). The 
normalisations in the original survey results are corrected using equation (47) and the total luminosity function (Table 2). 


# 

Survey 

- 5 loglO h 


±(e) 

^0 

(e) 

(10-2 Mpc-°) 

it' 

(10® h Lq Mpc-°) 

1 

SSRS2 

normalisation corrected 

-19.63+°;“ (bj)“ 

-1.00 ±0.09 

0.30 ± 0.08 

0.44 ± 0.08 

0.64 ±0.19 

0.41 ±0.09 

0.60 ±0.19 

2 

2dFGRS 

normalisation corrected 

-19.61 ±0.09 (bj) 

-0.74 ±0.11 

0.34 ± 0.06 

0.90 ±0.09 

0.82 ±0.17 

0.75 ±0.10 

0.68 ±0.15 

3 

2dFGRS 

normalisation corrected 

-19.59 ±0.05 (bj) 

-0.54 ±0.02 

0.40 ± 0.04 

0.99 ±0.05 

1.01 ±0.13 

0.79 ± 0.05 

0.81 ±0.11 


“ The B{0) magnitude of SSRS2 was converted to bj magnitude using bj = B — 0.28(i? — V) (Blair & Gilmore 1982; Norberg et al. 
2001) and assuming B{0) = B (Alonso et al. 1993) and mean (B — V) = 0.94 (Norberg et al. 2001). 


sizes) used for spectroscopic observations in various galaxy 
redshift surveys are significant sources of systematic errors 
in galaxy classifications. More generally, those techniques 
will give correct results only up to the accuracy of the as¬ 
sumed relation between the morphological types and (the 
parameters used to describe) the spectra of galaxies. In¬ 
deed, applying spectrum-based classification techniques to a 
sample of 4000 2dFGRS galaxies whose morphological types 
are visually known, Madgwick (2003) hnds that the rates of 
successfully classifying early-type and late-type galaxies by 
spectrum-based classification techniques range from ~ 60%- 
80%. Nonetheless, the results by Folkes et al. (1999) and 
Madgwick et al. (2001) are some of the few relatively more 
reliable results that are available. The values of the Schechter 
parameters for the local early-type galaxy population ob¬ 
tained by Folkes et al. (1999) and Madgwick et al. (2001) 
can be found in Table 3. In Table 3, the corrected early-type 
normalisations using equation (47) can also be found next 
to each survey result. While the corrected normalization is 
higher for the SSRS2 result, it is essentially not changed ei¬ 
ther for the Folkes et al. (1999) result or the Madgwick et 
al. (2001) result. This implies that the total luminosity den¬ 
sities measured by Folkes et al. (1999) and Madgwick et al. 
(2001) are consistent with our estimate while that measured 
by M98 is lower than our estimate. In section 4, we will use 
for our analyses the corrected SSRS2 (M98) early-type LF 
and the early-type LF directly measured by Madgwick et al. 
(2001). 

Since the total luminosity density is dominated by late- 
type galaxies, we assume that the faint-end slope and the 
characteristic absolute magnitude for late-type galaxies are 
the same as those for all galaxies, namely we take = 
-1.22 ± 0.02 and All’ll = -19.69 ± 0.04 (bj). Under this 
prescription, the normalisation for the late-type population 
is the same as the ratio of the luminosity density of the 
late-type population to the total luminosity density. Hence 
the inferred normalisations for the late-type population are 

= 1.20 ±0.11 10~^ Mpc“® for the M98 type-specihc 
LFs, = 1.13±0.09 Mpc“® for the Folkes et al. (1999) 
type-specihc LFs, and = 1.03±0.08 10“^ Mpc“® for 
the Madgwick et al. (2001) type-specihc LFs. 

Finally, it is required that galaxy evolution (of relevance 
to lensing) be incorporated in the analyses of statistical 
lensing because the lensing probability is sensitive to the 
redshift of the potential lensing galaxy (Turner, Ostriker & 
Gott 1984). Since the lensing rate is dominated by early-type 


galaxies, the most relevant issue of galaxy evolution in sta¬ 
tistical lensing is whether there is any signihcant change in 
the comoving number density of early-type galaxies at most 
likely lens redshifts (i.e. z ~ 0.6) compared with the present 
epoch as would be the case in certain rapid-evolution models 
(e.g. Baugh, Gole, & Frenk 1996). There exist in the litera¬ 
ture two qualitatively diherent kinds of observational results 
on the formation and evolution of early-type galaxies, each 
of which supports an extreme view of the formation and 
evolution of early-type galaxies; namely, an early formation 
by the monolithic collapse of protogalactic gases followed by 
passive evolution (Eggen, Lynden-Bell & Sandage 1962) or a 
continuous formation through merging of subunits in a hier¬ 
archical structure-formation model (Larson 1974; White & 
Rees 1978). Galaxy count results that are consistent with the 
early formation and passive evolution of early-type galaxies 
appear to be favoured (Im et al. 1996, 1999, 2002; Schade et 
al. 1999; Totani & Yoshii 1998; Lilly et al. 1995) while some 
results would support a rapid evolution of early-type galax¬ 
ies since an intermediate to high redshift (Fried et al. 2001; 
Kauffmann, Chariot & White 1996). Results from the anal¬ 
ysis by the fundamental plane of a sample of gravitational- 
lens early-type galaxies (Kochanek et al. 2000) would not be 
consistent with a rapid evolution of early-type galaxies but 
would be consistent with the early formation and passive 
evolution hypothesis. A further review of the current status 
on the formation and evolution of early-type galaxies can 
be found in Peebles (2002), who argues for the early forma¬ 
tion and passive evolution of early-type galaxies based on 
numerous lines of observational studies such as the colour- 
magnitude relation and the fundamental plane. Following 
Peebles (2002) and the authors mentioned above who ar¬ 
gue for passive evolution of early-type galaxies, it will be 
assumed in this study that the local early-type LF has not 
evolved in the sense that the characteristic comoving num¬ 
ber density and the faint-end slope of the early-type LF are 
unchanged in the lookback times. 


4 RESULTS 

The results on cosmological parameters and global proper¬ 
ties of galaxies that we derive in this section are under the 
following two key assumptions: 

• A galaxy as a gravitational lens is assumed to be well- 
approximated by a singular isothermal ellipsoid mass model. 
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Furthermore, dynamical normalisations of lenses are only 
considered for the oblate and the prolate cases. 

• The early-type comoving characteristic number density 
and faint-end slope are assumed to be unchanged from 2 ~ 1 
to the present epoch. 

We emphasize that the above assumptions are valid at least 
as a hrst-order approximation according to the majority of 
current evidence (section 2.1; section 3.3). The parameters of 
our statistical lensing model are summarized in section 2.1.3. 
Our goals are to constrain cosmological parameters and the 
early-type characteristic velocity dispersion (and the late- 
type characteristic velocity dispersion) and the mean pro¬ 
jected mass ellipticity of (early-type) galaxies from the ob¬ 
served properties of statistical lensing (given in section 3.1). 
To meet our goals, we need to fix the rest of the parame¬ 
ters because there exist degeneracies in the total parameter 
space with the observed properties of statistical lensing. Re¬ 
cent advanced observations have determined the fixed pa¬ 
rameters with relatively small uncertainties, as given in sec¬ 
tion 3.2, section 3.3, and below. Nonetheless, uncertainties 
in some parameters remain possible sources of systematic er¬ 
rors in the derived results and are discussed in sections 4.3 
and 5.3. 

The fixed parameters are as follows: (1) both the early- 
type and the late-type LF faint-end slopes and characteristic 
number densities as given in section 3.3 (note, however, that 
two alternative choices of the type-specific LFs are used; see 
below); (2) the Faber-Jackson exponent ypj = 4.0 (Bernard! 
et al. 2001) and the Tully-Fisher exponent 7 tf = 2.9 (Tully 
& Pierce 2000) (these parameters affect relatively weakly 
constraints on galaxy velocity dispersions and consequently 
cosmological parameters within their present uncertainties); 
(3) the parameters describing the redshift distribution and 
the differential number-flux-density relation for the CLASS 
sources as given in section 3.2 (systematic uncertainties aris¬ 
ing from the uncertainties in these parameters will be esti¬ 
mated separately). As simplihcations of our analyses, we use 
a common mean projected mass ellipticity e and a common 
relative frequency of the oblatenesses "Pobi both for the early- 
type and the late-type populations. The mean ellipticity and 
the intrinsic shape distribution of late-type galaxies have lit¬ 
tle effect on our derived results (other than the dependence 
of the derived late-type characteristic velocity dispersion on 
these parameters), so we simplify our analyses by using the 
above prescriptions. The parameters that are allowed to be 
free are as follows: (1) the cosmological parameters; matter 
density fJm, density of a classical cosmological constant Ua 
or dark energy density Qx with its constant equation of state 
w; (2) the early-type and the late-type characteristic velocity 
dispersions and (3) the mean projected mass el¬ 
lipticity of (early-type) galaxies e and the relative frequency 
of the oblatenesses "Pobi in the intrinsic shape distribution 
of (early-type) galaxies. We emphasize that it is not only 
justified but also very sensible to treat the characteristic ve¬ 
locity dispersions (in particular as free parameters and 
determine them from the observed image separations in the 
final statistical sample. This is because the observed radio¬ 
image separations are very accurate and as can be seen from 
equation (35), image separation has a strong sensitivity on 
velocity dispersion but little on cosmological parameters in 
particular within the ranges to be considered below. 


While we allow the early-type and the late-type charac¬ 
teristic velocity dispersions to be free, a characteristic cen¬ 
tral stellar velocity dispersion (at the present epoch) for 
the early-type population can be obtained through a mea¬ 
sured Faber-Jackson relation from the early-type charac¬ 
teristic absolute magnitude, and a characteristic maximum 
rotation speed 0 present epoch) for the late- 

type population can also be obtained through a measured 
Tully-Fisher relation from the late-type characteristic abso¬ 
lute magnitude. The early-type characteristic central stel¬ 
lar velocity dispersion and the late-type characteristic max¬ 
imum rotation speed (obtained in the above ways) can then 
be used to estimate the early-type and the late-type char¬ 
acteristic velocity dispersions, although the conversions are 
model-dependent. Below we use these methods to obtain 
independent estimates of the early-type and the late-type 
characteristic velocity dispersions. Such estimates can be 
compared with the derived values from the statistical lensing 
model fitting. Furthermore, we will use the estimated late- 
type characteristic velocity dispersion as a prior constraint 
for a model fitting in section 4.1. This choice of model fitting 
is motivated because there is only one conhrmed late-type 
lensing galaxy in the final CLASS statistical sample that 
can be used to constrain the late-type characteristic veloc¬ 
ity dispersion. 

The local early-type characteristic absolute magnitude 
of — 51ogj^Q h = —19.61 along with B = bj + 0.26 

(section 3.3) corresponds to a characteristic central stellar 

(e) _1 

velocity dispersion of = 192 ± 34 km s through 

—Mb = a + 9(logjQ Uc — T3) with a mean A — 19.5 ± 0.7 
for both ellipticals and SO galaxies (expressed for h = 1; 
de Vaucouleurs & Olson 1982). By htting the observed light 
distributions and the observed velocity dispersion prohles 
of 37 elliptical galaxies using singular isothermal mass dis¬ 
tributions, Kochanek (1994) finds that the central stellar 
velocity dispersion of an elliptical galaxy is a good mea¬ 
sure of the velocity dispersion of a singular isothermal mass 
model. Following Kochanek (1994), we thus estimate that 
afl « '^c*^o ~ 192 ± 34 km s“^. Alternatively, using the 
‘Faber-Jackson’ relation for the dark matter velocity disper¬ 
sion (i.e., the uum-L relation) shown in Fig. 4 of Kochanek 
(1994), we find that the above early-type characteristic ab¬ 
solute magnitude corresponds to an early-type characteristic 
velocity dispersion of = 185^® km s“^. 

We adopt the following Tully-Fisher relation to obtain a 
characteristic velocity dispersion afg for the late-type pop¬ 
ulation (Tully & Pierce 2000): 

- = 7.27[logio(2nWx) - 2.5] + 20.11, (48) 

where Umax is the maximum rotation speed (i.e. for the edge- 
on view) and is the extinction-corrected absolute B 
magnitude (i.e., the magnitude for the face-on view), which 
is given by 

M^ = M^‘’-{1.57 + 2.75[logio(2n«x)-2.5]}log(a/fe),(49) 

where is the uncorrected absolute magnitude and a/b 
is the apparent major-to-minor axis ratio. From M^^g(bj) — 
51ogjQ h = —19.69 ± 0.04 along with B — bj + 0.26 (sec¬ 
tion 3.3) and h = 0.72 ± 0.08 (Freedman et al. 2001), 
we obtain the characteristic maximum rotation speed of 
km s ^ using equations (48) and (49) and 
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taking a mean a/b — 2. If we assume an SIS mass distri¬ 
bution for the halo of a late-type galaxy, this characteristic 
maximum rotation speed corresponds to the characteristic 
velocity dispersion of = 134^^^ g km s“^ 

As discussed at length in section 3.3, while recent obser¬ 
vational determinations of the total galaxy luminosity func¬ 
tion give a converging result, the decomposition of the total 
luminosity function by morphological type is less certain. For 
this uncertainty, we do our analyses for the two choices of 
the type-specific galaxy LFs, i.e. the normalisation-corrected 
SSRS2 type-specific LFs and the 2dFGRS type-specific LFs 
(see Table 3 and section 3.3). The SSRS2 type-specific LFs^^ 
will be our standard choice because their galaxies were clas¬ 
sified through visual identifications (see section 3.3), and the 
2dFGRS type-specific LFs will be our alternative choice. No¬ 
tice that the SSRS2 early-type LF has a significantly lower 
characteristic number density and a significantly steeper 
faint-end slope compared with the 2dFGRS early-type LF 
[i.e. nl';^/(10”^ Mpc“®) = 0.64±0.19 vs. 0.99±0.05 and 
a = -i.OO ± 0.09 vs. -0.54 ± 0.02]. 

4.1 Constraints on cosmological parameters 

We first consider cosmological models with a classical cos¬ 
mological constant in which the cosmological parameters to 
be constrained are the matter density flm and the cosmolog¬ 
ical constant density Ha. In other words, we only consider 
a fixed constant equation of state for dark energy, w = — 1 . 
Notice, however, that even if we allow w along with flm and 
Qa to be varied simultaneously, the constraints on the flm- 
Ha plane would change little because the best-fitting value 
of w is —1 (with the prior constraint u; > — 1 imposed) 
at least for cosmological models that are near the flat line 
rim + riA = 1. We then consider flat cosmological models 
with a dark energy whose equation of state is assumed to be 
constant and satisfies the energy condition w > — 1. In this 
case, the cosmological parameters to be constrained are the 
matter density nm(= 1 — fla:) and the constant dark-energy 
equation of state. Notice here that because of degeneracy on 
the rim-flA plane along a line that is nearly orthogonal to 
the flat line for a given value of w, it is useful to impose the 
constraint flm + Qa = 1 . 

Fig. 6 shows confidence limits in the parameter plane 
spanned by flm and Qa- At each grid point on the plane the 
value of the (equation 44) is determined by minimizing 
the function over the nuisance parameters (i.e., all the other 
free parameters). The four contours respectively correspond 
to 68 %, 90%, 95%, and 99% confidence limits for one param¬ 
eter, namely Ax^ = 1.0, 2.71, 4.0, and 6.63 from the global 
minimum value of the Xmin- Notice that the adopted 
method for estimating confidence limits on model parame¬ 
ters is the so-called (profile) likelihood ratio method which is 
one of the widely practiced methods [see, e.g., Lupton (1993) 
and references therein]. This method gives more conservative 
(i.e. larger) confidence intervals than the frequently prac¬ 
ticed method in which nuisance parameters are held fixed 
at their global maximum-likelihood values. The confidence 
intervals obtained by the likelihood ratio method also closely 

Hereafter, ‘SSRS2 type-specific LFs’ always refer to the 
normalisation-corrected SSRS2 type-specific LFs. 


resemble realistic confidence intervals obtained by project¬ 
ing the n-dimensional confidence region onto the parameter 
subspace of interest. However, while computation of the like¬ 
lihood distribution in the n-dimensional (in the present case 
5-dimensional) parameter space is costly even for a moder¬ 
ately grained grid, distribution of likelihood in one or two- 
dimensional parameter subspace can be computed easily us¬ 
ing an efficient optimization routine such as the downhill 
simplex method (see Press et al. 1992). Here and through¬ 
out in this paper all the confidence limits on model param¬ 
eters are based on the likelihood ratio method. However, 
we also compute a likelihood distribution in a pseudo full¬ 
dimensional grid. By the likelihood distribution in a pseudo 
full-dimensional grid we shall mean the likelihood distri¬ 
bution in the subspace spanned only by the most crucial 
parameters where each likelihood value is obtained by fix¬ 
ing the rest of the parameters at their maximum-likelihood 
values. For the case under consideration, the most crucial 
three parameters are flm, I2 a, and ai'^\ The pseudo full¬ 
dimensional likelihood distribution is available at the Jodrell 
Bank Observatory web page for gravitational lenses and the 
CL ASS. The web page also includes (profile) likelihood 
distributions in parameter planes. 

Fig. 6 (a) and (b) are the results based on the SSRS2 
type-specific LFs while Fig. 6 (c) and (d) are those based on 
the 2dFGRS type-specific LFs. For the results of Fig. 6 (b) 
and (d), the late-type characteristic velocity dispersion is 
constrained by the value obtained from the Tully-Fisher re¬ 
lation (see above). 

Fig. 7 shows the ui-flm plane in flat cosmology. As in 
Fig. 6 , four different results are presented in Fig. 7 depending 
on the type-specific LFs and the assumption on the late-type 
characteristic velocity dispersion. In Table 4, we summarize 
the 68 % and the 95% confidence limits on the cosmological 
parameters. For the models considered in Fig. 6 and Fig. 7, 
the fitted values of the early-type characteristic velocity dis¬ 
persion are somewhat different for the different early-type 
luminosity functions. Gonfidence limits on galactic parame¬ 
ters including the early-type characteristic velocity disper¬ 
sion are presented next in section 4.2. 


In deriving constraints on galactic parameters, cosmological 
parameters are varied to minimize the (equation 44) at 
each grid point. Thus, the derived constraints on galactic 
parameters depend on the allowed ranges of cosmological 
parameters. In this section, we only consider flat cosmology 
(i.e. flm -l-fla: = 1) with 0 < flm < 1 and fix the dark-energy 
equation of state w at w = —1 (which is the best-fit value 
when it is allowed to vary with the constraint w > — 1 ). 

Fig. 8 shows confidence limits in the parameter plane 
spanned by the early-type characteristic velocity dispersion 
and the late-type characteristic velocity dispersion . 

The six different results in Fig. 8 are for the three different 
cases of the intrinsic shape distribution of galaxies, i.e. all 
oblate, one half oblate and the other half prolate, and all 
prolate, each for the SSRS2 and the 2dFGRS type-specific 

http://www.jb.man.ac.uk/research/gravlens/numerical/numerical.html 


4.2 Constraints on global properties of galaxy 
populations 
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(a) SSRSS/No prior (b) SSRS2/Gaussian prior on 



Figure 6. Likelihood regions in the Qm-fiA plane. The four contours, respectively, correspond to the 68%, 90%, 95%, and 99% confidence 
limits using the Ax^ (or, likelihood ratio) statistic for one parameter. The four illustrated models are as follows. Models (a) and (b) are 
based on the SSRS2 type-specific LFs while models (c) and (d) are based on the 2dFGRS type-specific LFs. For models (a) and (c), 
no prior constraints are imposed on the nuisance parameters while for models (b) and (d), = 134)j)jQ g km s“^ (see section 4) in 

conjunction with the dynamical normalisation for the equal frequencies of oblates and prolates is used as a prior constraint assuming 


LFs. Fig. 9 shows the parameter plane spanned by the mean 
projected mass ellipticity of galaxies e and the early-type 
characteristic velocity dispersion for the same six cases 
considered in Fig. 8. In Table 5, we summarize the 68% and 
the 95% confidence limits on the galactic parameters. 


4.3 Systematic effects of varying the parameters 
for the source population 

The results presented in sections 4.1 and 4.2 were derived 
assuming the present best estimates of the redshift distri¬ 
bution and the differential number-flux-density relation for 
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(a) SSRS 2 /N 0 prior (b) SSRS2/Gaussian prior on u, 



(=Px/Px) ^ ^=Px/Px) 


(c) 2dFGRS/No prior (d) 2dFGRS/Gaussian prior on 



'W {=PjPx) ^ i=Px/Px^ 


Figure 7. Likelihood regions in the m-Qm plane with the prior that the universe is flat for the same models illustrated in Fig. 6. 


Table 4. Constraints on cosmological parameters for the four models illustrated in Fig. 6 and Fig. 7. 


Model 

Ha - 1.2nm 
MLEi 

68% limit 

95% limit 

nm(= 1 - Ha) 
MLEi 

68% limit 

95% limit 

w 

MLEi 

68% limit 

(a) 

0.40 

(-0.16,0.73) 

(-1.17,0.92) 

0.31 

(0.17,0.58) 

(0.08,1.06) 

-1 

< - 0.55 

(b) 

0.12 

(-0.50,0.50) 

(-1.52,0.74) 

0.40 

(0.24,0.68) 

(0.13,1.14) 

-1 

< -0.41 

(c) 

0.40 

(-0.10,0.71) 

(-0.98,0.90) 

0.31 

(0.17,0.55) 

(0.09,0.96) 

-1 

< -0.58 

(d) 

0.17 

(-0.36,0.52) 

(-1.21,0.75) 

0.38 

(0.23,0.62) 

(0.13,1.00) 

-1 

< -0.49 


^ Maximum likelihood estimate which is the value minimizing the given by equation (44). 
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(1) SSRS2/Half Oblate, Half Prolate (2) 2dFGRS/Half Oblate. Half Prolate 



140 160 180 200 220 240 260 140 160 180 200 220 240 260 

(km (km s~^) 


(3) SSRS2/A11 Oblate (4) 2dFGRS/All Oblate 



140 160 180 200 220 240 260 140 160 180 200 220 240 260 

(km s~^) (km s~^) 


(5) SSRS2/A11 Prolate (6) 2dFGRS/All Prolate 



140 160 180 200 220 240 260 140 160 180 200 220 240 260 

(km s~^) (km s~^) 


Figure 8. Likelihood regions in the plane with the prior that the universe is Hat, w = —1, and 0 < f!m < 1- Models (1), (3), 

and (5) are based on the SSRS2 type-specific LFs while models (2), (4), and (6) and based on the 2dFGRS type-specific LFs. For models 
(1) and (2), one half of the galaxies of each type are assumed to be oblate while the other half are assumed to be prolate. For models 
(3) and (4) all galaxies are assumed to be oblate while for models (5) and (6) all galaxies are assumed to be prolate. 
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Table 5. Constraints on galactic parameters in flat cosmology for the six models illustrated in Fig. 8 and Fig. 9. 


Model 

(km s ^ 
MLEi 

) 

68% limit 

95% limit 

(km s ^ 
MLEi 

) 

68% limit 

95% limit 

e 

MLEi 

68% limit 

95% limit 

(1) 

198 

(180,220) 

(162,248) 

117 

(101,136) 

(86,160) 

0.42 

(0.28,—) 

(0.17,-) 

(2) 

181 

(166,199) 

551 , 220 ) 

119 

5o3,139) 

57,163) 

0.42 

5 . 28 ,—) 

(0.17,-) 

(3) 

175 

(160,205) 

546,233) 

103 

^9,126) 

57 , 149 ) 

— 

5 . 34 ,—) 

(0.21,-) 

(4) 

160 

(147,185) 

536,207) 

105 

5i,128) 

58,152) 

— 

5 . 34 ,—) 

(0.21,-) 

(5) 

208 

(190,232) 

573,262) 

123 

5o6,144) 

5i,169) 

0.35 

5 . 25 , 0 . 51 ) 

(0.14,-) 

(6) 

191 

(175,209) 

562,233) 

126 

5o8,146) 

53 , 172 ) 

0.35 

5 . 25 , 0 . 51 ) 

(0.14,-) 


^ Maximum likelihood estimate which is the value minimizing the given by equation (44). 



z rjg 

Figure 10. Systematic effects of varying (a) the mean source redshift, z and (b) the \dN/dS\ slope for 5 < 30 mjy, rj 2 , on the maximum 
likelihood value of Qm in flat cosmology with ui = — 1. The thin vertical solid lines, respectively, correspond to the present best estimates 
of z and 172 while the thin vertical dashed lines correspond to their present uncertainties. 


the CLASS sources as given in section 3.2. In other words, 
the present uncertainties in the redshift distribution and 
the differential number-ffux-density relation for the CLASS 
sources were not incorporated in the results in sections 4.1 
and 4.2. Here we quantify the effects of varying the param¬ 
eters used to describe the CLASS source population. The 
source parameters that can have significant impacts on the 
derived results on the cosmological parameters are the mean 
redshift, 0 and the differential number-flux-density relation 
slope for flux densities lower than 30 mJy, r ]2 (section 3.2). 
Figure 10 (a) and (b) show respectively the effects of varying 


the parameters « and 772 on the maximum likelihood value 
of flm in Hat cosmology with w = — 1. 


5 DISCUSSION 

5.1 Understanding the sensitivities of the model 
parameters to the observed properties of 
statistical leasing 

The likelihood function [equation 41; or the function 
(equation 44)] for statistical lensing is a complicated func- 
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tion of the parameters of the statistical lensing model (sec¬ 
tion 2.1.3). Moreover, the data are characterized by many 
properties (section 3) which control the likelihood function. 
For these reasons, it is not straightforward to see how the 
cosmological and galactic parameters are constrained by the 
data. Here we discuss the sensitivities of the model parame¬ 
ters to the properties of the data including numerical tests. 

To help us see how the model parameters are con¬ 
strained by the data, we rewrite the function (equa¬ 
tion 44) as 

= Xu + Xl, (50) 

where Xu Xl s-r® respectively the contributions from the 
unlensed sources and the lensed sources, which are given by 


JVu 


2 

Xu 


-El- 




and 




2 

XL 


-El-) 




(51) 


(52) 


where the integrated lensing probability Pm"l(fc) and the 
differential lensing probability are given by equa¬ 

tions (42) and (43). From equations (51) and (52), we see 
that increasing (decreasing) lensing probabilities increases 
(decreases) Xv but decreases (increases) Xl-^^ The proce¬ 
dure of fitting the model parameters to the data is to min¬ 
imize the sum of the oppositely behaving Xu Xl by 
varying the lensing probabilities. 

The lensing probabilities depend on many factors in¬ 
cluding galactic and cosmological factors and the magnifi¬ 
cation bias factor which depends on the source properties 
(section 2.1.2). In other words, the lensing probabilities are 
determined by the combination of all those factors. Accord¬ 
ingly, the lensing probabilities can be varied in a number 
of different ways of varying the factors. For example, the 
lensing probabilities can be decreased by reducing the cos¬ 
mological factors (« D{0, z)‘^\d£/dz\), by decreasing a*, or 
even by reducing the cosmological factors and increasing ct* 
(in the last case, one of course has to reduce the cosmologi¬ 
cal factors sufficiently enough to compensate for increasing 
(7*). However, different combinations of the factors are dis¬ 
tinguished by the differences in their capabilities of fitting 
the various properties of the data. For example, combining 
a very large ct* and a very low value of the cosmological fac¬ 
tors can give a ‘correct’ lensing rate, but such a combination 
is discarded (i.e. it gives a relatively larger x^ value) because 
it does not fit the observed image separation distribution. 

Roughly speaking, for a given dynamical normalisation 
(section 2.1.1), the observed image separations constrain the 
characteristic velocity dispersions and ai"'* and the ra¬ 
tio of the numbers of doubles and of quadruples constrains 
the surface density ellipticity. For these constrained galac¬ 
tic factors and a magnification bias factor determined from 
the estimated source properties, the cosmological factors are 
constrained mostly by the ratio of the numbers of the lensed 


Quast & Helbig (1999) used example figures to demonstrate 
this point [see their Figs. (Al) and (A2)]. 


sources and of the unlensed sources and (much) less signif¬ 
icantly by the redshifts of the lenses and of the sources. In 
practice, the galactic and cosmological factors are, of course, 
simultaneously varied to minimize the total x^- 


In Table 6, we present the calculated values of the Xu 
(equation 51) and Xl (equation 52) along with the fitted key 
parameters and the predicted lensing rate for the four differ¬ 
ent models shown in Fig. 7 and Table 4 and two additional 
models to be compared. For models (a), (b), and (b'), the 
SSRS2 type-specific LFs are used while for models (c), (d), 
and (d'), the 2dFGRS type-specific LFs are used. Notice that 
for models (b) and (d) the late-type characteristic velocity 
dispersion is constrained using a Gaussian prior, whose con¬ 
tribution to the total x^ iSj however, not shown in Table 6. 
For the two additional models (b') and (d'), the character¬ 
istic velocity dispersions are fixed at = 225 km s“^ and 
= 145 km s“^, which are the preferred values in several 
previous analyses of statistical lensing (e.g. Kochanek 1996a; 
Falco, Kochanek, & Munoz 1998; Quast & Helbig 1999; Hel¬ 
big et al. 1999). As can be seen in Table 6, galaxy popula¬ 
tions following the 2dFGRS type-specific LFs with the fitted 
characteristic velocity dispersions [Model (c)] have the low¬ 
est value for the x^- Compared with Model (c), Model (a) in 
which galaxy populations follow the SSRS2 type-specific LFs 
have a larger value for the x^ by Ax^ = 1.53. The larger x^ 
value for Model (a) compared with Model (c) is mostly due 
to the increase in the value of Xl- The implication for this 
is that Model (a) does not fit the observed image separation 
distribution as well as Model (c) does. Likewise, other mod¬ 
els presented in Table 6 have larger x^ values mainly because 
of their larger values of the Xl compared with Model (c). In 
particular, Models (b') and (d'), in which the characteristic 
velocity dispersions are held fixed, have significantly larger 
Xl values compared with Model (c) because the predicted 
mean image separations are discrepant with the observed 
mean image separation for those models. To show the above, 
we present in Fig. 11 the predicted image separation distri¬ 
butions for Models (a), (b'), (c), and (d'). In Fig. 11 are 
also displayed the model predicted mean image separations 
and the observed mean image separation for the multiply- 
imaged systems whose image separations are used in the fit 
(see section 3.1). Finally, notice that all the models in Ta¬ 
ble 6 predict correct lensing rates despite the differences in 
their capabilities of fitting the observed image separation 
distribution. This is, of course, the consequence of adjusting 
the value of to fit the observed lensing rate for the fitted 
(or given) galactic parameters in each model. 


From the above examination of the example models, we 
can summarize how the galactic and cosmological parame¬ 
ters are constrained by the data as follows. For the given 
type-specific LFs, the observed image separations constrain 
the characteristic velocity dispersions through the Xl term 
and then the observed lensing rate constrains cosmological 
parameters through the sum of the Xl s-nd the Xu terms. 
Since the lensing rate can be fitted by adjusting cosmologi¬ 
cal parameters for any given galactic parameters, the lensing 
rate alone cannot constrain cosmological parameters unless 
one has a priori accurate knowledge of the galactic param¬ 
eters. 
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Table 6. Relative values and predicted lensing rates for various models of statistical lensing. The x^ is given by x^ = Xu 
Xu (equation 51) and Xl (equation 52) are respectively the contributions from the unlensed sources and the lensed sources. Model (a), 
(b), (c) and (d) are the four models illustrated in Fig. 7 and Table 4, here for the case of flat cosmology with a classical cosmological 
constant. Notice that for the above models, the early-type and the late-type characteristic velocity dispersions were fitted. On the other 
hand, for Model (b^) and (d^) the velocity dispersions were fixed at = 225 km s“^ and = 145 km s“^ (see section 5.1). As for 
Model (a) and (b), the SSRS2 type-specific LFs were used for Model (b') while as for Model (c) and (d), the 2dFGRS type-specific LFs 
were used for Model (d'). 


Model 

r2m 

(= 1-flA) 

(km s“^) 

(km s ^) 

Xu 

xl 

X^ = Xu + Xl 

Ax^ 

pa 

(a) 

0.31 

197.8 

116.7 

25.17 

186.96 

212.12 

1.53 

1/711 

(b) 

0.40 

204.8 

130.2 

25.46 

187.14 

212.60 

2.01 

1/703 

(b') 

0.65 

225.0 

145.0 

25.89 

188.41 

214.30 

3.71 

1/691 

(c) 

0.31 

180.8 

119.0 

25.08 

185.52 

210.59 

= 0 

1/713 

(d) 

0.38 

185.5 

130.8 

25.12 

185.83 

210.95 

0.36 

1/712 

(d') 

0.89 

225.0 

145.0 

25.71 

189.70 

215.41 

4.82 

1/696 


^ The observed lensing rate is 1/689. 
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Figure 11. Image separation probability distributions for A0 > 0.3 arcsec of multiply-imaged CLASS sources at redshift Zs = 2 with 
lens redshifts 0.3 ^ zi < 1 (assuming that the unmeasured lens redshifts are also in this range). Theoretical distributions predicted by 
four models which can be found in Table 6 are overplotted on the observed distribution (i.e. the histogram) of the 9 multiply-imaged 
sources in the final CLASS statistical sample whose image splittings are known (or likely) to be caused by single galaxies. The small arrow 
marks the observed mean image separation, namely (A0) class = 1.13 arcsec, while the four vertical lines denote the predicted mean 
image separations of the four models. Notice that the mean image separations predicted by Model (a) and (c), in which the early-type 
and the late-type characteristic velocity dispersions were fitted, agree well with the observed mean image separation (see section 5.1). The 
mean image separations predicted by Model (b') and (d'), in which the early-type and the late-type characteristic velocity dispersions 
were fixed at = 225 km s“^ and = 145 km (see section 5.1), are discrepant with the observed mean image separation. The 
distribution predicted by Model (c), in which the 2dFCRS type-specific LFs were used, agrees best with the histogram (this explains 
why the model has the lowest x^ value), whose bins, however, consist only of small numbers of data points. 
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5.2 Comparison with previous results in 
statistical leasing 

Previous analyses of statistical leasing were based on 
an optically-selected sample, partial samples from radio- 
selected CLASS sources and other radio-selected sources, 
and galaxy luminosity functions derived prior to the 2dF- 
GRS and the SDSS observations. The optically-selected 
sample used for previous analyses is summarized in 
Kochanek (1996a) and consists of a combined sample of 862 
highly luminous QSOs, out of which 5 QSOs are multiply- 
imaged due to galactic mass-scale lenses. Radio-selected 
multiply-imaged sources used for previous analyses include 
4 multiply-imaged flat-spectrum radio sources (e.g. Helbig 
et al. 1999) out of a total of 6 JVAS multiply-imaged sources 
(Browne et al. 2003) and several other multiply-imaged 
sources including CLASS sources. 

Based on the optically-selected sample by Kochanek 
(1996a) and radio-selected samples, some previous analy¬ 
ses of statistical lensing (Kochanek 1996a; Falco, Kochanek, 
& Munoz 1998; Quast & Helbig 1999) obtained limits on 
Qa in flat cosmology (or, on Ha — Hm) that would rule 
out a high Ha- Most notably, Kochanek (1996a) obtained 
Ha ^ 0.66 a 95% confidence with a maximum likelihood 
value of Ha ~ 0. This result is significantly different from 
the results on cosmological parameters that are derived in 
section 4.1. To understand the results by Kochanek (1996a) 
and the others mentioned above, it is important to notice 
that (1) they used for the early-type characteristic veloc¬ 
ity dispersion = 225 ± 22.5 km s“^ from Kochanek 
(1994) as a prior information and (2) this prior value of 
was also consistent with the observed image separations of 
the multiply-imaged sources used by them for their adopted 
early-type luminosity function (Kochanek 1996a). The latter 
point is the consequence of combining two things. First, the 
mean image separation for the 5 optically-selected lens sys¬ 
tems is (A0)opticai = 1-60 arcsec which is significantly larger 
than the mean separation of (A6) class = 1-13 arcsec for 
the 9 multiply-imaged sources (or, (A0 )qlass ~ 
sec excluding the confirmed spiral lens system 0218-1-357) 
in the final CLASS statistical sample (section 3.1) whose 
image separations are (likely to be) due to single galactic 
potentials.Second, the fitted value of is a function of 
the early-type faint-end slope (for a fixed Faber-Jackson 
exponent 7fj) as seen in section 4.2, and the early-type LF 
used by the above authors has a relatively steep faint-end 
slope of = —1.00 ± 0.15 similar to that for the SSRS2 
early-type LF (section 3.3) and thus requires a relatively 
larger value of (compared with a shallower faint-end 
slope as in the 2dFGRS early-type LF). In passing, we men¬ 
tion that the local early-type characteristic number density 
estimated by Kochanek (1996a) agrees well with that 
for the SSRS2 early-type LF. 

Like the above authors, Helbig et al. (1999) adopted 
= 225 ± 22.5 km s“^ in deriving limits on cosmological 
parameters based on their JVAS sample. Helbig et al. (1999) 

The radio-selected multiply-imaged sources collected by 
Kochanek (1996a) and used to constrain have a mean sep¬ 
aration of (A6},.,jUio(Kochanek) ~ 1-6 which is nearly equal to 
{A0)optical* 


obtained —0.85 < Ha < 0.84 at 95% confidence with a max¬ 
imum likelihood value of Ha > 0 for a flat cosmology. This 
result is somewhat different from that by Kochanek (1996a) 
in the sense that it marginally favors a positive Ha. However, 
it is also different from the results of this paper (section 4.1) 
which favor a HA-dominated universe. Given the difference 
between the value of adopted by Helbig et al. (1999) 
and the values fitted to GLASS lenses (section 4.2), it is 
perhaps not surprising to find that the Helbig et al. (1999) 
result and our results are different. Yet, it would be of in¬ 
terest to assess and understand the difference between the 
Helbig et al. (1999) result and our results more concretely 
and clearly, especially since the JVAS sample used by Helbig 
et al. (1999) is a GLASS subsample just as the final GLASS 
statistical sample is. What if was fitted to the lenses 
in the Helbig et al. (1999) sample rather than being fixed 
at = 225 km s“^? What factors other than the value 
of (in particular, the more realistic statistical lensing 
model, the adopted profile likelihood ratio method of esti¬ 
mating parameters and their errors, and sample variance) 
might be attributed to the difference between the results of 
Helbig et al. (1999) and of this paper? First of all, we recall 
that there is no significant difference in galactic parame¬ 
ters other than between Helbig et al. (1999) and this 
paper. To help us see separately the effects of varying statis¬ 
tical lensing model, parameter and error estimation method, 
and sample, we consider Fig. 12 which shows four variants 
of Fig. 6 (a) (i.e., results based on the SSRS2 type-specific 
LFs). Fig. 12 (a) and (b) are the results based on the same 
data used by Helbig et al. (1999)^® but using our model. 
Fig. 12 (a) is the result obtained considering only early- 
type galaxies and taking = 225 km s“^ [i.e., adopting 
the same approach used by Helbig et al. (1999)]. To obtain 
Fig. 12 (a), we also assumed that oblates and prelates have 
equal frequencies (i.e. dynamical normalisation A ~ 1) and 
took a mean projected galaxy mass ellipticity e = 0.42 which 
gives the best-fit to the final CLASS statistical sample (see 
Fig. 9). Fig. 12 (b) is the result obtained by applying our 
new approach adopted in this paper to the Helbig et al. 
(1999) data. This means that both early-type and late-type 
populations were considered and their characteristic veloc¬ 
ity dispersions were fitted to the observed image separa¬ 
tions of the Helbig et al. (1999) 4 multiply-imaged systems. 
Fig. 12 (c) and (d) are the results based on two subsam¬ 
ples of the final CLASS statistical sample (section 3.1) and 
using the same model and approach adopted in this paper. 
Fig. 12 (c) is based only on the JVAS subsample within 
the final CLASS statistical sample which comprises 1749 
sources including three multiply-imaged sources 0218-1-357, 
1422-1-231, and 2114+022. Fig. 12 (d) is based on a subsam¬ 
ple of 2308 sources including four multiply-imaged sources 
0445+123, 0712+472, 2045+265, and 2319+051 which were 
randomly selected from the subsample of 7209 sources in- 


A noticeable difference in sample selection criteria between the 
Helbig et al. (1999) sample and the final CLASS statistical sample 
is that the lower limit of the fainter-to-brighter image flux-density 
ratio for doubly-imaged sources is 0.05 for the Helbig et al. (1999) 
sample. This was taken into account in deriving the results shown 
in Fig. 12 (a) and (b). 
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Figure 12. Likelihood regions in the plane for the SSRS2 type-specific LFs as in Fig. 6(a) based on previous and present 

subsamples from the entire JVAS and CLASS sources, shown to illustrate the effects of different statistical lensing models, different 
parameter and error estimate methods, and different samples. The four contours in each panel, respectively, correspond to the 68%, 
90%, 95%, and 99% confidence limits for one parameter, (a) & (b) are based on the Helbig et al. (1999) sample; (a) is the result of only 
fitting the lensing rate fixing early-type characteristic velocity at = 225 km s“^ whereas (b) is the result of simultaneously fitting 
both the lensing rate and the image separations. Both (c) and (d) are the results of simultaneously fitting both the lensing rate and the 
image separations based on subsamples of the final CLASS statistical sample: (c) is the result based on the subsample that comprises 
only 1749 JVAS sources whereas (d) is based on a random subsample that comprises 2308 non-JVAS sources. The comparison of (a) 
and Fig. 1 of Helbig et al. (1999) shows the effects of the different statistical lensing models. The comparison of (a) and (b) shows the 
effects of the different parameter and error estimate methods. The comparison among (b), (c), and (d) shows differences or similarities 
among samples. Finally, the comparison of (d) [or (c)] and Fig. 6(a) shows the effects of different sample sizes. See section 5.2 for a full 
discussion of the results. 
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eluding 10 multiply-images sources within the final CLASS 
statistical sample that are not JVAS sources. 

Let us briefiy discuss the implications of the results dis¬ 
played in Fig. 12. From Fig. 12 (a), we have Ha = 0.48lg 3g 
(68%) and O.dSlg'H (95%) for a fiat cosmology. This re¬ 
sult favors a positive Ha universe and the difference from 
the Helbig et al. (1999) result is not significant. Thus, we 
are led to conclude that the differences in the details of the 
statistical lensing model (in particular, cross sections and 
magnification biases) between this paper and Helbig et al. 
(1999) make little difference in deriving limits on cosmologi¬ 
cal parameters based on presently available data (i.e. within 
presently reachable accuracies). From Fig. 12 (b), we have 
Ha = 0 . 32 ^g' 7 | (68%) for a fiat cosmology. The somewhat 
lower maximum-likelihood estimate of Ha from Fig. 12 (b) 
compared with that from Fig. 12 (a) is the consequence of 
the early-type characteristic velocity dispersion larger than 
225 km s“^ fitted to the image separations that are due 
to early-type galaxies in the Helbig et al. (1999) sample 
(the mean image separation for 0414-1-054, 1030-1-074, and 
1422-1-231 is 1.64 arcsec). Another main difference between 
Fig. 12 (a) and (b) is that the confidence interval at a given 
confidence level is much smaller in Fig. 12 (a) although both 
were based on the same data. This is, of course, due to the 
fact that was fixed for Fig. 12 (a) while it was varied for 
Fig. 12 (b) based on the profile likelihood ratio method (see 
section 4.1). The above comparison of Fig. 12 (a) and (b) 
shows that different approaches of treating galactic param¬ 
eters can make significant differences in the derived values 
of cosmological parameters for the same data. The results 
displayed in Fig. 12 (c) and (d) are based on the profile likeli¬ 
hood ratio method as in Fig. 12 (b). We have Ha = 0.68lg'4g 
and0.64l°;22 at 68% confidence respectively from Fig. 12 (c) 
and (d). These values of Ha are very similar to each other 
but different from the value obtained from Fig. 12 (b). For 
example, the confidence interval of Ha implied by the Hel¬ 
big et al. (1999) data excludes at 68% confidence both of the 
maximum-likelihood estimates of Ha shown on Fig. 12 (c) 
and (d). In other words, even if not at a high statistical sig¬ 
nificance, the Helbig et al. (1999) data are inconsistent with 
two subsets of the the final CLASS statistical sample when 
all three of the data sets are analysed and interpreted in 
the same manner. This could imply that the Helbig et al. 
(1999) sample is incomplete within their specified sample 
definitions or is not as observationally well-defined as the fi¬ 
nal CLASS statistical sample. Alternatively, the Helbig et al. 
(1999) sample simply (as a rare occurrence) happens to in¬ 
clude the lensed system 0414-1-054 that has an image separa¬ 
tion of 2.09 arcsec and is a main factor behind the low value 
of Ha. Comparing Fig. 12 (d) [or Fig. 12 (c)] with Fig. 6 (a), 
we find that as the sample size becomes four times larger, 
the 68% confidence interval approximately halves consistent 
with the expectation for a Poisson error. 

In summary, Kochanek (1996a), Helbig et al. (1999) and 
the other authors above either found or adopted « 225 
km s“^, and for this and other adopted galactic parameters 
the observed lensing rates in the optical sample of QSOs and 
the JVAS samples of fiat-spectrum radio sources were most 
consistent with relatively low values of Ha. In particular, 
some of these authors produced results that even appeared 
to rule out a high value of Ha such as Ha = 0.7. 

Some previous analyses of statistical lensing were based 


on choices of galactic parameters that were different from 
those used by Kochanek (1996a), Helbig et al. (1999), and 
the others above. Cooray (1999) used choices of and 

that were similar to those used by Kochanek (1996a) 
but separated the early-type population into an elliptical 
(E) population and an SO population with a number ratio 
of E: SO = 1: 1.6 and estimated the characteristic velocity 
dispersions from his chosen LFs and Faber-Jackson relation 
as (T* = 2IOI45 km s“^ for the E population and ct, = 
194ltig km s“^ for the SO population. The Cooray (1999) 
estimates of the E and SO characteristic velocity dispersions 
correspond to an effective characteristic velocity dispersion 
for the entire early-type population that is > la lower than 
the estimate by Kochanek (1994). Not surprisingly, Cooray 
(1999) obtained a relatively high upper limit of Ha <0.79 
at 95% confidence in flat cosmology from the lensing rate in 
an incomplete subsample of ~ 6500 CLASS sources. Notice, 
however, that Cooray (1999) did not use the observed image 
separations in his analysis to test whether the chosen LFs 
and characteristic velocity dispersions were consistent with 
them or not. 

Choices of galactic parameters that are significantly dif¬ 
ferent from those used by Kochanek (1996a), Cooray (1999), 
and the others above were considered by Chiba & Yoshii 
(1999) and Cheng & Krauss (2001). In particular, based on 
the same lens systems used by Kochanek (1996a), Chiba & 
Yoshii (1999) found Ha ^0 at 98% confidence with a max¬ 
imum likelihood value of Ha ~ 0.7 because of their radi¬ 
cally different choices of galactic parameters. The galactic 
parameters preferred by Chiba & Yoshii (1999) and Cheng 
& Krauss (2001) can be characterized by a very shallow or 
inverted (i.e. positive) faint-end slope for the early-type LF. 
Furthermore, Cheng & Krauss (2001) argued for a shallow 
slope of 7 « 2 to 3 in the aom-L relation for the early-type 
galaxies, and Chiba & Yoshii (1999) used an early-type char¬ 
acteristic number density [i.e. an early-type LF by Loveday 
et al. (1992)[ that is a factor of « 2 to 3 lower than those 
in the SSRS2 and the 2dFGRS early-type LFs (Table 3). 
The shallow/inverted early-type faint-end slopes necessarily 
require lower values for ai^^^ for the given same image sepa¬ 
rations compared with a steeper slope (as noted above). The 
high value for Ha obtained by Chiba & Yoshii (1999) was the 
combined effect of the significantly lower values for the early- 
type characteristic velocity dispersion and number density 
than those used by Kochanek (1996a). Although the limits 
on Ha obtained by Chiba & Yoshii (1999) agree well with 
the results on cosmological parameters based on the final 
CLASS statistical sample (section 4.1), we note that (1) the 
early-type LF used by Chiba & Yoshii (1999) is discrepant 
with present more reliable early-type LFs (section 3.3) and 
(2) the statistical lensing samples used by Chiba & Yoshii 
(1999) (i.e. those collected by Kochanek 1996a) were defined 
using selection functions that are observationally not as re¬ 
liable as those used to define the final CLASS statistical 
sample and thus are more likely to be subject to unknown 
biases (which may explain the difference in the mean image 
separations between the final CLASS statistical sample and 
other samples). 

Compared with the observational information that was 
available to the previous analyses of statistical lensing (in 
particular those that put limits on cosmological parameters), 
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the presently available observational information used in this 
work is more abundant and more reliable. The final CLASS 
statistical sample is defined so as to satisfy the selection cri¬ 
teria (section 3.1) that are both observationally well-defined 
and reliable. This means that its biases are well understood 
and can be accurately taken into account in the analysis of 
the sample. [For comparison, optically-selected samples suf¬ 
fer from dust extinction which cannot be accurately taken 
into account for analyses of statistical leasing (see Falco et 
al. 1999).] The final CLASS statistical sample is the most 
significant piece of information that has been available only 
to this work. The total galaxy LFs derived from the 2dF- 
GRS and the SDSS observations and more recently derived 
type-specific LFs have now significantly reduced uncertain¬ 
ties in the early-type LF. The fact that the well-defined final 
CLASS statistical sample and the relatively more reliable 
early-type LFs were used and all the factors in the statis¬ 
tical leasing model were carefully taken into account based 
on relatively reliable data and with thorough understand¬ 
ing, makes the results on cosmological parameters reported 
in Chae et al. (2002) and here (section 4.1) much more reli¬ 
able than those from previous analyses. 

So far in this section, we have focused on cosmolog¬ 
ical parameters in the comparison with previous analyses 
of statistical lensing. We now compare our derived results 
on galactic parameters with those in previous analyses. As 
emphasized above, the mean image separation in the fi¬ 
nal CLASS statistical sample is significantly smaller than 
those in the optical QSO sample and the heterogeneously 
defined radio-selected samples (Kochanek 1996a). We recall 
that this was the main reason for the difference in the cos¬ 
mological constraints between Kochanek (1996a) and this 
work (section 4.1). The implication for the early-type char¬ 
acteristic velocity dispersion is that it is ~ 27 km s“^ lower 
than the estimate by Kochanek (1994) for the steep SSRS2 
early-type LF and ~ 44 km s“^ lower for the shallow 2dF- 
GRS LF.^® We point out that this fitted value of the early- 
type characteristic velocity dispersion is also consistent with 
our estimate of the early-type characteristic magnitude (sec¬ 
tion 3.3). Namely, our estimate of the early-type character¬ 
istic magnitude Mfg{B) is ~ 0.5 mag fainter than the one 
adopted by Kochanek (1996a) and it implies ~ 190 
km s“^ either by a Faber-Jackson relation (de Vaucouleurs 
& Olson 1982) or a (Jdm-L relation (Kochanek 1994) (see 
section 4). 

Our statistical lensing model incorporates the (differ¬ 
ential) lensing probabilities of specific image multiplicities 
for the multiply-imaged sources using the SIE lens model 
(rather than a circular lens). This allows us to constrain 
the mean projected mass ellipticity of galaxies from the 
multiply-imaged sources of well-defined image multiplicities. 
In section 3.1, we pointed out that 10 systems out of the 
13 lens systems in the CLASS statistical sample have well- 
defined image multiplicities for which single galactic poten¬ 
tials are (likely to be) only responsible. The 10 systems con¬ 
tain 6 doubly-imaged ones and 4 quadruply-imaged ones. 
The relative frequency of the doubles in the final CLASS sta¬ 
tistical sample is similar to (or slightly higher than) those 

This is based on the dynamical normalisation for the equal 
numbers of the oblates and the prelates, i.e. A(/) ^ 1. 


in the samples used by King & Browne (1996), Kochanek 
(1996b), and Keeton, Kochanek, & Seljak (1997). However, 
the relative frequency for the doubles is higher than that 
in the CLASS subsample used by Rusin & Tegmark (2001), 
in which there are 5 doubles and 7 quadruples. This dif¬ 
ference is for two reasons. First, in the time since Rusin 
& Tegmark (2001) did their analysis based on CLASS lens 
systems available at that time, the CLASS group has iden¬ 
tified two new lens systems 0445-1-123 (Argo et al. 2003) 
and 0631-1-519, both of which are doubly-imaged systems 
and are now included in the final CLASS statistical sample. 
Second, the observational criteria used to define the final 
CLASS statistical sample (section 3.1) are not only more 
stringent but more reliable than those used to define the 
statistical sample used by Rusin & Tegmark (2001). Based 
on their statistical sample and the SIE lens model, Rusin 
& Tegmark (2001) concluded that the predicted mean mass 
ellipticity of the early-type galaxies (tmass ~ 0.6) was much 
higher than an observed mean light ellipticity . In particular, 
Rusin & Tegmark (2001) found that the observed ellipticity 
distribution of the early-type galaxies in the Coma cluster 
was inconsistent with the relative frequencies of the dou¬ 
bles and the quadruples in their sample at 98% statistical 
significance. 

Based on the final CLASS statistical sample and the 
SIE lens model, we find that the best-fit mean mass ellip¬ 
ticity for the (early-type) galaxies^® is Cmass = 0.42 (Fig. 9; 
Table 5) and the 68% lower limit is 0.28 for the dynami¬ 
cal normalisation of the equal frequencies of the oblates and 
the prelates. For the dynamical normalisation of all oblates, 
the 68% lower limit is increased to 0.34. For the dynami¬ 
cal normalisation of all prelates, tmass = 0.357Q ^g at 68%. 
These results are consistent with the observed mean light 
ellipticity eiight,DBEP/GSS ~ 0.38 (Fig. 18 in Im et al. 2002) 
for the field early-type galaxies of the Deep Extragalactic 
Evolutionary Probe (DEEP) Gross Strip Survey (CSS) or 
Hight.Coma ~ 0.32 for the combined E and SO galaxies in the 
Coma cluster [Fig. 3(a) in Jprgensen & Franx 1994]. There¬ 
fore, the potential ‘ellipticity crisis’ for gravitational lenses 
that were noticed and investigated in the previous works 
(Rusin & Tegmark 2001; Kochanek 1996a; King & Browne 
1996) based on samples that are not as well-defined as the 
final CLASS statistical sample, has now significantly weak¬ 
ened. 


Rusin & Tegmark (2001), as in this work, did not use com¬ 
pound lenses in which multiple potentials are responsible for the 
image multiplicities. 

Recall that a common mean mass ellipticity for both early- 
type and late-type populations of galaxies is considered and fitted 
to the image multiplicities that are generated by both early-type 
and late-type (appropriately isolated) lensing galaxies. Since most 
of the lensing galaxies are early-type, the common mean mass 
ellipticity can be regarded as the mean mass ellipticity of early- 
type galaxies with little error. 
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5.3 Possible sources of systematic errors 

5.3.1 Mass profiles, intrinsic shapes, and velocity 
dispersions of galaxies 

In this work, we adopted the singular isothermal mass pro- 
hle and the axisymmetric intrinsic shapes (i.e., oblate and 
prolate) to model galaxy lenses. We also assumed that the 
distribution function of a galaxy depends only on the rela¬ 
tive energy and the angular momentum component parallel 
to the symmetry axis of the galaxy. With the adopted mass 
model and the assumption on the distribution function, we 
calculated, as functions of apparent axial ratio and line-of- 
sight velocity dispersion, the multiple-imaging cross sections 
averaged over all the combinations of intrinsic shapes and 
inclination angles that are consistent with the apparent ax¬ 
ial ratio for the assumed relative frequencies of oblates and 
prolates. The apparent axial ratio and the characteristic ve¬ 
locity dispersions (corresponding to the characteristic ab¬ 
solute magnitudes) were then constrained by the observed 
relative frequencies of image multiplicities and the observed 
image separations in the hnal CLASS statistical sample. We 
found that the values of the early-type characteristic veloc¬ 
ity dispersion obtained in the way described above agreed 
well with the values based on the early-type characteristic 
absolute magnitude and a Faber-Jackson relation or a udm- 
L relation (recall here that a singular isothermal mass profile 
was also assumed; see section 4). The values of the late-type 
characteristic velocity dispersion obtained from the statisti¬ 
cal lensing analyses were « Icr lower than the value based on 
the late-type characteristic absolute magnitude and a Tully- 
Fisher relation. 

How would the derived results from the statistical lens¬ 
ing analyses be affected if non-isothermal mass profiles, non- 
axisymmetric intrinsic shapes, and/or general three-integral 
distribution functions of galaxies were considered? Con¬ 
straints on cosmological parameters are not likely to be af¬ 
fected by non-isothermal mass proHles, non-axisymmetric 
intrinsic shapes, or general distribution functions according 
to the following arguments. For the given apparent axial ra¬ 
tio and line-of-sight velocity dispersion, the cross sections 
averaged over the allowed ensemble of models obviously 
depend on the assumed galaxy model. Hence, the derived 
characteristic velocity dispersions from the observed image 
separations depend on the mass profile, the allowed range 
of intrinsic shapes, and the assumption on the distribution 
function of the galaxy. (Dependences on the intrinsic shapes 
for the singular isothermal mass profile were illustrated in 
section 4.) However, while the derived characteristic veloc¬ 
ity dispersions are affected as the galaxy model is varied, 
the multiple-imaging cross sections are not as long as the 
projected surface density of the inner cylindrical region of 
the galaxy is similar to the isothermal surface density. In 
other words, the observed image separations will determine 
the multiple-imaging cross sections regardless of the details 
of the galaxy model as long as the image separation for the 
given multiple-imaging cross section is not signihcantly de¬ 
pendent on the image magnification ratio, similarly to the 
case for the singular isothermal profile in which the image 
separation is fixed by the cross section and independent of 
the magnihcation ratio. 


5.3.2 Local early-type LF: the characteristic number 
density and faint-end slope 

As discussed in section 3.3, there is no consensus among ob- 
servationally derived results for the early-type LF. We have 
used two choices of the early-type LF in this work, namely 
the steep SSRS2 choice and the shallow 2dFGRS choice. 
As seen in section 4, while the two different early-type LFs 
require different values of to fit the same image separa¬ 
tions in the CLASS statistical sample, the difference in the 
values of for the two early-type LFs is such that the two 
early-type LFs give similar results on cosmological param¬ 
eters (Fig. 6 & 7). However, the fact that the two different 
early-type LFs give similar results on cosmological parame¬ 
ters does not preclude the possibility that the presently de¬ 
rived constraints on cosmological parameters based on the 
two early-type LFs suffer from systematic errors due to the 
present uncertainties in the early-type LF. As pointed out in 
section 3.3, the main difficulty in the observational deriva¬ 
tion of the early-type LF lies in classifying large numbers of 
galaxies by morphological types. Systematically misclassify- 
ing galaxies will not only result in an error in the number 
density of early-type galaxies but also an error in the faint- 
end slope in the early-type LF. Given that the derived cos¬ 
mological parameters depend on both the early-type charac¬ 
teristic number density and faint-end slope, it is understood 
that the presently derived constraints on cosmological pa¬ 
rameters are potentially susceptible to unquantified system¬ 
atic errors due to the present uncertainties in the early-type 
LF. Nonetheless, since the total galaxy LF has been reliably 
determined, as long as the partition of the total LF into the 
type-specihc LFs is not drastically different from those con¬ 
sidered in this work, the conclusions drawn on cosmological 
parameters (section 4.1) would not be significantly affected. 


5.3.3 Evolution of early-type galaxies 

One key assumption made in this work is no evolution of 
early-type galaxies since 2 ~ 1 in the sense that the early- 
type characteristic comoving number density and faint-end 
slope are unchanged from 2 ~ 1 to the present epoch. In 
section 3.3, we pointed out that the hypothesis of the early 
(2 > 2) formation and passive evolution of early-type galax¬ 
ies is supported by many lines of observational studies. 

However, the evolution of early-type galaxies remains a 
relatively larger source of uncertainty for analyses of statis¬ 
tical lensing. First, the observational results that provided 
the most direct pieces of evidence for little or no evolution of 
the early-type characteristic comoving number density and 
faint-end slope were based on relatively small samples of 
galaxies and were less than abundant. Furthermore, lim¬ 
its on galaxy evolution derived from galaxy counts are de¬ 
pendent on the assumed cosmological model and thus, in 
a strict sense, should not be used for constraining cosmo¬ 
logical parameters via statistical lensing (see, e.g., Keeton 
2002). Second, there exist observational results in the lit¬ 
erature that have been used to claim a rapid evolution of 
early-type galaxies in which there is a significant reduction 
in the comoving number density of early-type galaxies at in- 
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termediate redshifts (0.3 ^2^1) compared with the present 
epoch. 

If future observational studies established a rapid evo¬ 
lution in the population of early-type galaxies since 2 ~ 1 so 
that the comoving number density of early-type galaxies at 
an intermediate redshift was significantly lower than that at 
the present epoch, the constraints on cosmological param¬ 
eters derived in this work would have to be adjusted; the 
adjustment would increase (decrease) dark energy density 
(matter density). Reversing the problem, statistical lensing 
can be used in the future to constrain the formation and evo¬ 
lution of early-type galaxies assuming a certain cosmological 
model. From the results of this work, the qualitative infer¬ 
ence on the formation and evolution of early-type galaxies 
would be that a flat ‘concordance cosmology’ with « 0.3 
(Ha ~ 0.7) favors an early (2> 1.5) formation of early-type 
galaxies followed by passive evolution. A detailed quantita¬ 
tive analysis of statistical lensing focusing on constraining 
the formation and evolution of early-type galaxies assuming 
the concordance cosmology will be presented in a forthcom¬ 
ing publication (Mao & Chae 2003). 

5 . 3.4 Redshift distribution of CLASS sources 

As shown in Figure 10(a), the derived value of (in a flat 
universe) is sensitive to the mean redshift of CLASS sources. 
Since the Marlow et al. (2000) measurement has a relatively 
large uncertainty due to the relatively small number of mea¬ 
sured spectroscopic redshifts, it will be important to obtain 
secure redshifts for more CLASS sources in the future to 
reduce the uncertainty arising from the current uncertainty 
in the mean redshift of CLASS sources. 

Six multiply-imaged sources in the Huai CLASS statisti¬ 
cal sample (Table 1) lack measured redshifts. While the un¬ 
known redshifts would not change the derived results from 
analyses of statistical lensing as long as they are not too close 
to their corresponding lens redshifts,^® it would be impor¬ 
tant to measure them not just as a matter of completeness 
but to make sure that no surprising results are missed. 


6 FUTURE PROSPECTS OF STATISTICAL 
LENSING 

As demonstrated in detail in this work, statistical properties 
of strong gravitational lensing allow us to probe indepen¬ 
dently cosmological models and global properties of galaxies. 
In principle, a sample of cosmologically distant sources that 
is observationally well-dehned and large enough can alone 

However, it is interesting to note that those who found a rapid 
evolution of early-type galaxies classified galaxies through photo¬ 
metric/spectroscopic information (e.g. Fried et al. 2001) while 
those who found little or no evolution of early-type galaxies clas¬ 
sified galaxies through morphological appearances or light profiles 
(e.g. Im et al. 2002). 

The uncertain redshift of 2045+265 can be an issue here. The 
system has the highest lens redshift among the measured lens 
redshifts and a relatively large image separation. Furthermore, 
there was an early report that the source redshift could be Zs = 
1.28 although this interpretation was very uncertain (Fassnacht 
et al. 1999). 


pin down both cosmological parameters (e.g. flm, LI a, and lo) 
and global parameters of galaxies (e.g. , early-type veloc- 

(s) 

ity function slope, cr; , spiral-type velocity function slope, 
and tmass) including their evolutionary properties. Since de¬ 
termining cosmological parameters and global properties of 
galaxies is a main and fundamental part of cosmology, it is 
important to have various independent methods. Statistical 
lensing provides such a method. 

The Cosmic Lens All-Sky Survey is the largest system¬ 
atic strong gravitational lens survey at the time of this writ¬ 
ing. A second Cosmic Lens All-Sky Survey (CLASS2) is be¬ 
ing planned and would aim to discover ten times as many 
multiply-imaged radio-loud sources as the CLASS did. A 
statistical sample from a CLASS2 project would allow us to 
measure cosmological and galactic parameters with three 
times smaller statistical errors compared with the values 
obtained in this work assuming that the errors are Pois- 
sonian [this appears to be the case from the comparison of 
Fig. 12 (c) or (d) with Fig. 6 (a) (see section 5.2) and we 
will deal with the issue regarding errors in a future publica¬ 
tion through an extensive Monte Carlo simulation of statis¬ 
tical lensing (Chae et al., in preparation)]. For example, this 
means that statistical lensing could independently come up 
with a value of flm (or Ll^) with an error of ~ 0.05 (if a flat 
universe is assumed) and values of and with errors 
of ~ 7 km s“^. 

It is almost certain that statistical lensing will not 
enable us to measure cosmological parameters with the 
precision that is achievable by cosmic microwave back¬ 
ground anisotropy measurements or supernovae observa¬ 
tions through space-based instruments for the present and 
the coming decade [e.g., Wilkinson microwave anisotropy 
probe {WMAP: Spergel et al. 2003); Planck surveyor^^; su¬ 
pernova/acceleration probe (SNAP)^^]. Nevertheless, inde¬ 
pendent results on cosmological parameters from statisti¬ 
cal lensing are of interest and value because the importance 
of cosmological parameters warrants demanding tests. How¬ 
ever, perhaps the greatest value of statistical lensing for the 
coming years lies in its unique probe of the velocity functions 
of galaxies and galaxy evolution. In particular, the velocity 
dispersions implied by the image separations of strongly- 
lensed systems are different from those measured spectro¬ 
scopically from galaxy surveys such as the SDSS, and statis¬ 
tical lensing probes galaxy evolution based on mass-selected 
samples of galaxies whereas studies from galaxy surveys are 
based on brightness-selected samples. Thus, statistical lens¬ 
ing provides valuable and unique constraints on theories of 
the formation and evolution of galaxies. 


7 CONCLUSIONS 

From the analyses of statistical lensing based on the final 
CLASS statistical sample, the current estimate of luminos¬ 
ity functions of galaxies per morphological type, and the sin¬ 
gular isothermal ellipsoid (SIE) lens model, we have reached 
the following conclusions: 


http://astro.esa.int/SA-general/Projects/Planck 
http://snap.lbl.gov 
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(i) The multiple-imaging cross sections of galaxies im¬ 
plied by the observed image separations are such that the 
observed ‘local’ distributions of galaxies require a cosmolog¬ 
ical model dominated by dark energy to be compatible with 
the observed multiple-imaging rate, under the assumption 
that field early-type (i.e. E/SO) galaxies existed with the 
present-epoch abundance and masses at intermediate red- 
shifts 0.3 < « < 1. In particular, the confidence region in the 

plane is in good agreement with that from Type la 
supernovae observations (Riess et al. 1998; Perlmutter et al. 
1999). This independent result further supports a flat con¬ 
cordance cosmology with ~ 0.3 (Ha ~ 0.7). However, 
if the early-type galaxy population underwent a significant 
evolution from a ~ 1 to the present epoch as in a certain 
hierarchical galaxy formation picture, the statistical lensing 
data used in this work would imply an even higher dark 
energy density. 

(ii) Assuming that one half of early-type galaxies 
are oblates and the other half are prelates, the fitted 
early-type characteristic line-of-sight velocity dispersion is 

(0.3 < « < 1) = 198t jg km s“^ (68%) if the early-type LF 
faint-end slope follows that of the ‘steep’ SSRS2 early-type 
LF (i.e. = -1), or cri'’'(0.3<2< 1) = 18lll^ km s-^ 

(68%) if it follows that of the ‘shallow’ 2dFGRS early-type 
LF (i.e. = —0.54). These results are consistent with 

the current estimate of the early-type characteristic abso¬ 
lute magnitude in the B band and the a-L relation, but 
are significantly lower than the previous result based on 
the optically-selected lens sample and the SSRS2-like steep 
early-type faint-end slope. 

(iii) Assuming that one half of early-type galaxies are 
oblates and the other half are prelates, the best-fit value for 
the mean projected mass ellipticity of early-type galaxies is 
0.42 and the 68% lower limit is 0.28. The inferred mean pro¬ 
jected mass ellipticity agrees well with the measured mean 
light ellipticity of early-type galaxies in ‘fields’ (i.e., random 
environments) or in clusters. 
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